Submitted to ApJ, June 5, 2008 

Preprint typeset using I4TgX style emulateapj v. 10/09/06 



HOW DO DISKS SURVIVE MERGERS? 
Philip F. Hopkins 1 , Thomas J. Cox 1 - 2 , Joshua D. Younger 1 , & Lars Hernquist' 

Submitted to ApJ, June 5, 2008 



oo 

o 

O 

(N 
•*-> 
O 

O 

(N 

S 1 

Oh- 

6 



(N 
> 

m 

o 

oo 
o 



X 



ABSTRACT 

We develop a general physical model for how galactic disks survive and/or are destroyed in mergers and 
interactions. Based on simple dynamical arguments, we show that gas primarily loses angular momentum to 
internal torques in a merger, induced by the gravity of the secondary. Gas within some characteristic radius, 
determined by the efficiency of this angular momentum loss (itself a function of the orbital parameters, mass 
ratio, and gas fraction of the merging galaxies), will quickly lose angular momentum to the stars sharing the 
perturbed host disk, fall to the center and be consumed in a starburst. We use a similar analysis to determine 
where violent relaxation of the pre-merger stellar disks is efficient on final coalescence. Our model describes 
both the dissipational and dissipationless components of the merger, and allows us to predict, for a given 
arbitrary encounter, the stellar and gas content of the material that will survive (without significant angular 
momentum loss or violent relaxation) to re-form a disk in the merger remnant, versus being dissipationlessly 
violently relaxed or dissipationally losing angular momentum and forming a compact central starburst. We 
test these predictions with a large library of hydrodynamic merger simulations, and show that they agree well 
(with small scatter) with the properties of simulated merger remnants as a function of merger mass ratio, 
orbital parameters, and gas distributions, in simulations which span a wide range of parameter space in these 
properties as well as prescriptions for gas physics, stellar and AGN feedback, halo and initial disk structural 
properties, redshift, and galaxy masses. We show that, in an immediate (short-term) sense, the amount of 
stellar or gaseous disk that survives or re-forms following a given interaction can be understood purely in 
terms of simple, well-understood gravitational physics, independent of the details of the ISM gas physics or 
stellar and AGN feedback. This allows us to demonstrate and quantify how these physics are in fact important, 
in an indirect sense, to enable disks to survive mergers, by lowering star formation efficiencies in low mass 
systems (allowing them to retain large gas fractions) and distributing the gas to large radii. The efficiency of 
disk destruction in mergers is a strong function of gas content - our model allows us to explicitly predict and 
demonstrate how, in sufficiently gas rich mergers (with quite general orbital parameters), even 1 : 1 mass-ratio 
mergers can yield disk-dominated remnants, and more realistic 1:3-1:4 mass-ratio major mergers can yield 
systems with < 20% of their mass in bulges. We discuss a number of implications of this modeling for the 
abundance and morphology of bulges as a function of mass and redshift, and provide simple prescriptions for 
the implementation of our results in analytic or semi-analytic models of galaxy formation. 

Subject headings: galaxies: formation — galaxies: evolution — galaxies: active — galaxies: spiral — cosmol- 
ogy: theory 



1. INTRODUCTION 

In the now established "concordanc e" ACDM cosmol - 
ogy, structure grows hierarchically (e.g. White & Rees 1978), 
making mergers and interactions between galaxies an es- 
sential and inescapable process in galaxy formation. In- 
deed, mergers are widely believed to be responsible for 
the morphologie s of spheroids (bulges in disks and el- 
liptical galaxies; Toomrej Il977l) . and observations find re- 
cent merger remnants in considerable abundance in the 
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local univer se (Schweizer 1982; Lake & Dressier 1986; 
Dovon et alj|1994t IShier & Fischer! 1 1998b Uames et all 1999 : 



lMaUeretalJl200a iFakhouri & m3EQQ1 IStewart et al.| |2007 



and references therein) and observations (e.g.lLin et al|2004 ; 
iBarton et alJl2007l:IWoods et al.ll2006HWoods & Gellerf2007() 



Genzeletal. 2001; Tacconi et al 
2007; Rothberg &Joseph| |2004 



2002; iDasvra et all |200€ , 
120061) as well as e.g. 
faint shells and tid al features common a round apparently 
"normal" galaxies dMalin & Carter! 1 198011 19831; ISchweizeri 
Il980t ISchweizer&Seitzerl 119921: ISchweizeri 1 19961) . which 
are thought to be signatures of galaxy collisio ns (e.g. 

iHernquist & Ouinrill988t IHernquist & Spergelll 19921). 

From both theoretical grounds (Ostriker & Tremaine 1975; 
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it appears that "minor" mergers of mass ratios < 1 : 10 are 
ubiquitous (there are almost no galaxies without mergers of at 
least this mass ratio in the last few Gyr), and moreover a large 
fraction (~ 1/2) of the ~ galaxy population is observed 
and expected to have experi enced a "major" merger (mass ra- 
tio < 1 : 3) s ince z ~ 2 - 3 dLotz et alj|2008l; iBell et alj|2006t 
iBridge et all l2007t iLin et alj 120081: iKartaltepe et all 120071) . 
With increasing redshift, kinematic and morphological indica- 
tions of recent, violen t disturbance in disk-dominated galaxies 
appear more frequent (|Hammer et al j2005l ; lFlores et all2006t 
iPuech et al]l2007ll2008h . 

Far from there not being enough mergers to explain the 
abundance of bulges and ellipticals, this has led to the concern 
that there may be far too many mergers to explain the survival 
and abundance of galactic disks in the context of our presen t 
understanding of galaxy formation. Toomre & ToomreTd 19721) 
were among the first to point out that mergers are capable 
of dramatically altering the morphologies of disks, trans- 
forming them into elliptical galaxies. Although their ne- 
glect of the importance of dissipational star formation and 
gas dynamics in the mergers led to some controversy (e.g. 
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IOstrikerifT98Qt ICarlberdll986t IGudj3[1987I: lKormendvlfT989h . 
it is now increasingly well-established that major mergers be- 
tween spiral galaxies (similar to those observed locally and 
at z < 2-3) with gas fractions comparable to those observed 
yields remnants in good agreement with essentially all ob- 
served properties of low and intermediate-mass local ellip- 
tical galaxies (e.g. morphologies, shapes, sizes, kinematics, 
densities, colors, black hole properties, fundam ental scal- 
ing relations, stellar populations, and halo gas: | Her nquist 



1989; Barnes & Hernquist 1991, 1996; Hernquistet al. 1993; 
Mihos & Hernquistl Il994bl 119961; iDi Matteo et aTTboOsf 



Naab et all l2006t Uesseit et al.1 120071; ICox et all 12006a 



Robertson etail l2006bt ISpringel et al I |2005at iBurkert et ail 
20071; iHopkins et alj20073l 2b08a. g.b.clldT). 

Many intermediate and low-luminosity "cusp" ellipticals 
(encompassing ^ 80 - 90% of the mass density in el- 
lipticals) contain significant embedded dis ks (perhaps all 
such e llipticals, given pr ojection effects; see Ferrares e et alj 
Il994t lLauer et al.1 120051) . and they form a continuous se- 
quence with most SO galaxies, known to have promi 
nent stellar (and e ven gaseous) disks dKormendvl 1 19851 



1994; Lauer et al. 1995; Faber et al 



Ferrarese et al. 2006; Emselle m et al.l 



Bender etal.1 19921; iFerrarese et al. Il994t iKormendv et al 



997t lKormendylll999t 



2007). Indeed, the ex- 
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istence of embedded disks in simulated merger remnants is 
critical to matching the properties described above. 

A wide variety of observations including stellar pop- 
ulations and star formation histories (e.g. | Bender et al.l 
Il989t iTrager et all 120001; iMcDermid etail 120061) and kine- 
matic and structural analysis of recent merger remnants 
( Schweizer et al.ll983tlSchweizedl983tlSchweizer & S eitze 
I1992L I2007t iHibbard & Yunlll999[ iRoth berg & J osephll200 
demonstrate that most of these disks are not accreted in the 
standard cosmological fashion after the spheroid forms - they 
must somehow survive the merger or form very quickly there- 
after from gas already in and around the galaxies. There- 
fore, despite the destruction of a large portion of a stellar disk 
in major mergers, some disk must survive mergers, and the 
amount that does so is a critical component determining many 
of the photometric and kinematic properties of even bulge- 
dominated and elliptical galaxies. 

Moreover, "minor" mergers - at least those with mass ra- 
tios < 10 : 1 (below which the difference between "merger" 
and accretion becomes increasingly blurred) - are not gen- 
erally believed to entirely destroy disks, but they are almost 
an order of magnitude more frequent than major mergers and 
as such may pose a more severe a problem for disk survival. 
In the ACDM cosmology, and from observed satellite frac- 
tions, it is unlikely than any disk (let alone a large fraction 
of disk galaxies) with a significant stellar age has survived ~ 
5-10 Gyr without ex periencing a merger of mass ratio 10 : 1 
or larger. Simulations dOuinnll 984; Ouinn & Goodma n! 19861; 



Ouinn etalJ Il993t iHernq uist & Mihosl 119951; I Walker etail 
1996t IVela zquez & White 1999; Naab & Burked 12003b 



Bournaud et alJ2005H Younger et al.ll2007ll2008h and analytic 
arguments dOstriker & Tremaindl975l ; lToth & Ostrikerl 19921 ; 
Sellwo od et alJ Fl998) suggest that gas-poor minor mergers 
can convert a considerable fraction of a stellar disk into bulge 
and cause significant perturbation ("puffing up" via dynamical 
heating) to the disk. The observed coldness of galactic disks 
sugges ts that this may be a severe problem: iToth & Ostrikerl 
( 1992) argued that large disks such as that in the Milky Way 
could not have undergone a merger of mass ratio S, 10 '■ 1 in 
the last ~ lOGyr. More recently e.g. lStewart et al.ld2007l) and 



iHammer et al.l d2007l) emphasized that the tension between 
these constraints and the expectation in CDM models that a 
number of such mergers should occur implies either a deficit 
in our understanding of hierarchical disk formation or a chal- 
lenge to the concordance cosmological model. 

Given the successes of the ACDM model on large scales, 
and the increasing observational confirmation that disks do 
undergo (and therefore must somehow survive) a large num- 
ber of mergers, it is likely that the problem lies in our 
(still relatively poor) understanding of disk galaxy forma- 
tion. This has led to a great deal of focus on the problem 
of forming realistic disks in a cosmological context, with 
many different attempts and debate on the missing elements 
necessary to produce disks in simulations. Various groups 
have argued that self-consistent treatment of gas physics 
and star formation along with implementation of feedback 
of different kinds is necessary, along with greatly improved 



numerical resolution dWeil et aljfl9 98; Sommer- Larsen et al 
1999, 2003; Thacker & Couchman 2000, 2001; Abadi et al 



20031: iGovernato et al.l I200H 120071: iRo bertson eTail 12004 
Okamoto et al.ll2005l : IScannapieco et al.ll2008l) . in order to en- 
able disks to survive their expected violent merger histories 
without completely losing angular momentum and transform- 
ing into systems that are too compact and have too much bulge 
mass (relative to real observed disks) by z = 0. 

It has been kn own fo r some time (s ee e. 

IHernquist & Barnes! 1199 It iBarnes & Hernquistl fl99 
that (even without any feedback) some fraction of the gas 
in even a major merger of two disks can survive and form 
new, embedded disks in the remnant - i.e. despite the prob- 
lems outlined above, disks are not necessarily completely 
destroyed in mergers. However, early studies of this were 
restricted to cases with low gas content (/ gas < 10% in the 
progenitor disks), most of which was rapidly consumed in star 
formation, yielding small remnant dis ks in strongly bulge- 
domin ated remnants. In seminal w ork, [Springel & Hernquist 
(2005) and Rob ertson et al.1 ((2006a) showed that, in idealized 
merger simulations with significant stellar feedback to allow 
the stable evolution of extremely gas rich disks (/ gas ~ 1), 
even a major merger can produce a disk-dominated remnant. 
This ha s since been confirmed in fully cosmological simu- 
lations dGovernato et alj [2007). Together with other recent 
investigations (see references above), these works have led to 
the growing consensus that a combination of strong stellar 
feedback and large gas content is essential to the survival of 
disk galaxies. 

A large number of open questions remain, however. How, 
exactly, does feedback allow disks to survive mergers? What 
are the most important physics? Does it require fine-tuning of 
feedback prescriptions? How might things vary as a function 
of galaxy mass, redshift, gas content, merger orbits, and envi- 
ronment? Fundamentally, should this be expected for typical 
cosmological circumstances, or are these cases pathological? 

The ambiguity largely owes to the fact that there is no deep 
physical understanding of how disks survive or re-form after 
mergers and interactions. It has only just become possible to 
conduct simulations with the requisite large gas fractions, and 
thus far theoretical explanations have largely been restricted 
to phenomenological analysis, with continued efforts to im- 
prove resolution and sub-resolution prescriptions. Moreover, 
without a full model for how disks behave in interactions, 
these simulations cannot be placed into the broader context 
of the emergence of the entire Hubble sequence (for example 
asking the question, are the disks in lenticulars and embedded 
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disks in ellipticals survivors of their pre-merger disks? Are 
they re-accreted? What determines how large they are? What 
is the key physics that gives rise to realistic embedded disks, 
leading to bulge-dominated galaxies with kinematic and pho- 
tometric properties similar to those in the real universe?) or 
within a fully cosmological context. 

The resolution requirements for full models of disk for- 
mation are severe - limiting any attempt to properly simu- 
late a cosmological box and still achieve the resolution nec- 
essary to reliably model a disk population - and so mod- 
els of the population of disks, largely semi-analytic, are 
forced to adopt simplified and un-tested prescriptions for 
the behavior of disks in mergers. This, in turn, has led 
to other well-known problems in modeling disk populations 
(even where prescriptions can ensure no artificial angular 
momentum losses); even when the cumulative (morphology- 
independent) galaxy mass function is correctly predicted at 
the low-mass end, semi-analytic models widely overproduce 
the relative abundance of low-mass spheroids and underpro- 
duce disks (even when satellites, which have other associated 
model uncertainties, are removed from consideration; see e.g. 
Somerville et al.ll200lL 120081; ICroton et al.ll2006t iBower et all 
20061: Ide Lucia & Blaizotl 12007b . Lacking a proper, physi- 
cally motivated understanding of how low-mass or gas-rich 
disks may or may not survive mergers, attempts to address 
this problem in the models have been purely phe nomenolog- 
ical a nd involve arbitrary prescriptions (see e.g. Ko da et al] 
l2007h . 

Motivated by these concerns, in this paper we develop a 
physical, dynamical model for how disks survive and are de- 
stroyed in mergers and interactions. We show that, in an im- 
mediate (short-term) sense, the amount of stellar or gaseous 
disk that survives or re-forms following a given interaction 
can be understood purely in terms of simple, well-understood 
gravitational physics. Knowing these physics, we develop an 
analytic model that allows us to accurately predict how much 
of a given pre-merger stellar and cold gas disk will survive a 
merger, as a function of the merger mass ratio, orbital param- 
eters, pre-merger cold gas fraction, and mass distribution of 
the gas and stars. We compare these predictions to the results 
of a large library of hundreds of hydrodynamic simulations 
of galaxy mergers and interactions, spanning a wide param- 
eter space in these properties as well as prescriptions for gas 
physics, stellar and AGN feedback, halo and initial disk struc- 
tural properties, redshift, and absolute galaxy masses. Our 
numerical experiments confirm that the analytic scalings ac- 
curately describe the behavior and bulge formation/disk de- 
struction in mergers over the entire dynamic range surveyed, 
and confirm that the parameters not explicitly included in our 
model do not systematically affect either the mean predictions 
or the scatter of simulations about those predictions. This al- 
lows us to understand the mean behavior of systems with dif- 
ferent orbits and mass ratios, as well as why systems with 
large gas fractions can form little bulge in even major merg- 
ers. 

This is possible because gas, in mergers, primarily loses 
angular momentum to internal gravitational torques (from the 
stars in the same disk) owing to asymmetries in the galaxy 
induced by the merger. Hydrodynamic torques and the direct 
torquing of the secondary are second-order effects, and very 
inefficient. Once gas is drained of angular momentum, there 
is little alternative but for it to fall to the center of the galaxy 
and form stars, regardless of the details of the prescriptions 
for star formation and feedback (these may change things at 



the ~ 10-20% level by blowing out some of the gas, but 
they cannot fundamentally alter the fact that cold gas with no 
angular momentum will be largely unable to form any sort 
of disk, or the fact that a galaxy's worth of gas compressed 
to high densities and small radii will inevitably form a large 
mass in stars). But if the systems are sufficiently gas-rich, 
then there is little stellar material sharing the disk to torque on 
the gas in the interaction, and little or no angular momentum 
is lost. 

Feedback can dramatically alter the ability of a disk to sur- 
vive in a cosmological sense: by allowing galaxies to retain 
large gas fractions (as opposed to no-feedback scenarios, in 
which cold gas in a disk is usually quickly converted into 
stars), they are more gas-rich when they undergo interactions, 
allowing them to avoid angular momentum loss for the reason 
above. Moreover, we show that in detail (owing to the res- 
onant structure of interactions), it is really gas within a cer- 
tain radius of the stellar disk that is drained of angular mo- 
mentum. The commonly-invoked stellar wind feedback then 
enables cosmological disk survival in a second fashion: by 
redistributing gas out to large radii, it prevents angular mo- 
mentum loss and allows rapid re-formation of disks after a 
merger. Independent of any tuning, our model allows us to 
quantify the disks expected as a function of interactions of ar- 
bitrary properties, and to physically, explicitly quantify what 
the requirements are for feedback, in a cosmological scenario, 
to enable disk survival. 

In § |2] we describe our library of gas-rich merger simula- 
tions, which we use to test our physical model for disk de- 
struction and survival. In §[3] we demonstrate the existence of 
genuine disks in remnants of even major mergers and briefly 
consider their properties, and compare methods to separate 
the disks and bulges in merger remnants. In §[4] we consider 
the question of how these disks form in and survive merg- 
ers: we identify the key components of any merger remnant 
in § 14.11 highlighting that these disks originate from a combi- 
nation of undestroyed pre-merger stellar disks and ga s wh ich 
avoids angular momentum loss in the merger. In § 14.21 we 
discuss how, in detail, that angular momentum loss proceeds. 
We use this, in § 14.31 to build a physical model for how angu- 
lar momentum loss proceeds in mergers and predict the sur- 
viving disk content of merger remnants: we model and test 
how this depends on the gas content of the pre-merger disks 
(§ 14. 3. U . and the orbital parameters (§ 14.3.2b and mass ratio 
(§ 14. 3. 31 of the encounter. We generalize to first passage and 
fly-by encounters (§ I4.3.41 i and demonstrate that (for other- 
wise fixed conditions at the time of an encounter) our conclu- 
sions are purely dynamical, independent of feedback physics 
or details in our treatment of e.g. star formation and the ISM 
gas physics (§ 14.3.51 ). although we use our model to deter- 
mine exactly how these choices can have dramatic indirect 
consequences for disk survival (by altering the state of sys- 
tems leading into a merger). We discuss some exceptions and 
pathological cases in § 14.3.61 and relate our r esults t o the long- 
term secular evolution of barred systems in § 14.3.71 In §|3J we 
outline how these results can and should be applied in ana- 
lytic and semi-analytic models of galaxy formation, and give 
appropriate prescriptions derived from our numerical experi- 
ments. Finally, we summarize our results and discuss some 
of their cosmological implications and applications to other 
models and observations in § [6] 

Throughout, we assume a Hm = 0.3, ft\ = 0.7, Hq = 
70kms" 1 Mpc~ 1 cosmology, but this has little effect on our 
conclusions. 
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2. THE SIMULATIONS 

Our simulation s were perform ed with the parallel TreeSPH 
code GADGET-2 ( Springe ll2005l). employing the fully conser- 
vative formulation (Springel & Hernquist 2002) of smoothed 
particle hydrodynamics (SPH), which conserves energy and 
entropy simultaneous ly even when smoothing lengths evolv e 
adaptively (see e.g., iHernquistl 119931: lO'Shea et al] 120051) . 
Our simulations account for radiative cooling and incor- 
porate a sub-resolution model of a multiphase interstellar 
medium ( ISM) to describe star form ation and supernova 
feedback (Springel & Hernquisl l2003l) . Feedback from su- 
pernovae is captured in this sub-resolution model through 
an effective equation of state for star-forming gas, en- 
abling us to stably evolve disk s with arbitrary gas fractions 
(see, e.g. ISpringel et al.[l2005bl: ISpringel & Hernquistl 120051: 
iRobertson et alj|2006allcl) . This is described by the parameter 
leos, which ranges from q eos = for an isothermal gas with 
effective temperature of 10 4 K, to q eos = 1 for our full multi- 
phase model with an effective temperature ~ 10 5 K. We have 
also compared with a subset of simulati ons which adopt the 
star form ati on fe edback prescription from Mi hos & Hernquistl 
(ll994allbUcl 119961) . in which the ISM is treated as a single- 
phase isothermal medium and feedback energy is deposited 
as a kinetic impulse. We examine the effects of these choices 
in § 14.3.51 and find they are minimal. 

Likewise, although they make little difference to the anal- 
ysis here, supermassive black holes are usually included at 
the centers of both progenitor galaxies. These black holes 
are represented by "sink" particles that accrete gas at a rate 
M estimated from the local gas density and sound speed us- 
ing an Eddington-limited prescription based on Bondi-Hoyle- 
Lyttleton accretion theory. The bolometric luminosity of the 
black hole is taken to be Lboi = e r Mc 2 , where e r = 0.1 is the 
radiative efficiency. We assume that a small fraction (typi- 
cally rj 5%) of Lboi couples dynamically to the surrounding 
gas, and that this feedback is injected into the gas as thermal 
energy, weighted by the SPH smoothing kern el. This fraction 
is a fre e parameter, which we determine as in Di Matt eo et al] 
(2005) by matching the observed Mbh - o relation. For now, 
we do not resolve the small-scale dynamics of the gas in the 
immediate vicinity of the black hole, but assume that the time- 
averaged accretion rate can be estimated from the gas proper- 
ties on the scale of our spatial resolution (roughly 20 pc, in 
the best cases). While the black holes can be indirectly im- 
portant, owing to their feedback ejecting gas into the halo and 
thus preserving it from star formation until the final merger, 
we find that, for a given gas content at the time of the actual 
merger, our results are unchanged in a parallel suite of simu- 
lations without black holes. 

The proge n itor galaxy models are described in 
Springel et al. (2005b), and we review their properties 
here. For each simulation, we generate two stable, isolated 
disk galax i es, ea ch with an extended dark matter halo with a 
iHernq uist (1990) pr ofile, motivated by cosm ological simula- 
tions (iNavarro et alJl!996t iBusha et al]l2005l) . an exponential 
disk of gas and stars, and (optionally) a bulge. The galaxies 
have total masses M vil - = v4/(10G#[z]), with the baryonic 
disk having a mass fraction m& = 0.041, the bulge (when 
present) having = 0.0136, and the rest of the mass in dark 
matter. The dark mat ter halos are ass i gned a concentration 
parameter scaled as in Roberts on et al .1 (l2006d) appropriately 
for th e galaxy mass and redshift following iBullock et al j 
d2001l) . We have also varied the concentration in a subset of 



simulations, and find it has little effect on our conclusions 
(because the central regions of the galaxy are, in any case, 
baryon-dominated), insofar as they pertain to disk survival in 
mergers (it has been demonstrated that halo concentrations 
are important for e.g. the exact sizes and velocity scalings of 
disks, and our predicted disk sizes scale accordingly). The 
initial disk scale-length is computed based on an assumed 
spin parameter A = 0.033, chosen to b e near the mode in th e 
A distribution measured in simulations (IVitvitska et al.l l2002). 
and the scale-length of an initial bulge (when present) is set 
to 0.2 times this. 

Typically, each galaxy initially consists of 168000 dark 
matter halo particles, 8000 bulge particles (when present), 
40000 gas and 40000 stellar disk particles, and one BH parti- 
cle. We vary the numerical resolution, with many simulations 
using twice, and a subset up to 128 times, as many particles. 
We choose the initial seed mass of the black hole either in 
accord with the observed Mbh-c relation or to be sufficiently 
small that its presence will not have an immediate dynamical 
effect, but we have varied the seed mass to identify any sys- 
tematic dependences. Given the particle numbers employed, 
the dark matter, gas, and star particles are all of roughly equal 
mass, and central cusps in the dark matter and bulge are rea- 
sonably well resolved. 

We consider a series of several hundre d simulations of 
collidin g galax ies, desc r ibed in detail in Rob ertson et al j 
(2006b]i) and |Cox_etaLl (2006a b). We vary the numerical 
resolution, the orbit of the encounter (disk inclinations, peri- 
center separation), the masses and structural properties of the 
merging galaxies, presence or absence of bulges in the pro- 
genitor galaxies, initial gas fractions, halo concentrations, the 
parameters describing star formation and feedback from su- 
pernovae and black hole growth, and initial black hole masses. 

The progenitor galaxies have virial velocities V wu = 
55,80,113,160,226,320, and 500kms-\ and redshifts z = 
0,2, 3, and 6, and the simulations span a range in final 
spheroid mass Mbh ~ 10 8 - 10 13 M Q , covering essentially the 
entire range of the observations we consider at all redshifts, 
and allowing us to identify any systematic dependences in our 
models. We consider initial disk gas fractions (by mass) of 
/ gas = 0.05, 0.1, 0.2, 0.4, 0.6, 0.8, and 1.0 for several choices 
of virial velocities, redshifts, and ISM equations of state. The 
results described in this paper are based primarily on simula- 
tions of equal-mass mergers; however, we examine in § 14.3.31 
how our results scale with mass ratio in mixed encounters, 
down to mass ratios ~ 1 : 10 or so, below which (as we show) 
the encounters have little noticeable eff ect. In detail, t h e sim- 
ulations studied there are described in lYounger et al] (120081) 
and constitute a complete subset of permutations of our stan- 
dard galaxy models with mass ratios uniformly sampling the 
range 1 : 1 to 1 : 8. As in our larger set of 1 : 1 mergers, at 
each mass ratio we systematically survey the effects of dif- 
ferent absolute galaxy mass, orbital parameters, and disk gas 
fraction (resulting in a typical ^30-40 simulations spanning 
the full range of orbital parameters and gas fractions of inter- 
est, around each mass ratio ~ 1 : 1, 1:2, 1:4, and 1 : 8). We 
have considered more limited studies of minor mergers where 
we vary e.g. the ISM equation of state, redshift, initial disk 
structural properties; as we find in our studies of these param- 
eters in the larger suite of equal-mass mergers, they make no 
significant difference to our conclusions here. 

Once built, pairs of galaxies are placed on parab olic or- 
bits (motivated by cosmological simulations; see e.g. [Benson] 
2005; Khochfar & Burkert 2006) with the spin axis of each 
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TABLE 1 
Disk Orientations 
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disk specified by the angles 6 and <fi in standard spherical 
coordinates. Table Q] lists the orientations in different rep- 
resentative orbits we have s ampled . The particular choice 
of orbits follows ICox et a l. (2006b); there are seven ideal- 
ized mergers (cases b-h) that represent orientations often seen 
in the literature (for example, case h, where all the angu- 
lar momentum vectors of th e disks and o rbit are initially 
aligned), the rest (i-p) follow Barnes (1988) by selecting un- 
biased initial disk orientations according to the coordinates 
of two oppositely directed tetrahedrons. These orbits are 
id entical to those conside red in various other studies, such 
as iNaab & Burkertl d2003l). For our s eries of orbits of var- 
ious mass ratios from Younger et al.l (12008b (where, in mi- 
nor mergers, the inclination of the secondary is less important 
than that of the primary) we survey the inclination of the pri- 
mary in a systematic sense, considering all our mergers with 
0i = 0, 30, 90, 150, 180° (cas es m00 0-ml80). We examine 
the effect of orbits in detail in § 14.3.21 and find that for random 
orbits, the differences are quantifiable but not strong - patho- 
logical orbits (such as the aligned case h above) are discussed 
in § l4.3.6l (these pathological cases often, in fact, are the most 
efficient at destroying disks). We have also tested our pre- 
dicted scalings with limited subsets of simulations that vary 
the pericent ric passage distance and the energy of the orbi t, 
described in IRobertson et all (l2006al) and lCox etall (l2006bl) . 
and find that our estimates are robust to these variations. 

Each simulation is evolved until the merger is complete and 
the remnants are fully relaxed, typically ~ 1 — 2 Gyr after the 
final merger and coa lescence of the BH s. We then analyze the 
remnants following ICox et all d2006bl) . in a manner designed 
to mirror the methods typically used by observers. For each 
remnant we project the stars onto a plane as if observed from 
a particular direction (we consider 100 viewing angles to each 
remnant, which uniformly sample the unit sphere). When we 



plot projected quantities such as R e , a, and V c , we typically 
show just the median value for each simulation across all ~ 
100 viewing directions. The sightline-to-sightline variation 
in these quantities is typically smaller than the simulation-to- 
simulation scatter, but we explicitly note where it is large. 

3. THE EXISTENCE OF DISKS IN MAJOR MERGER REMNANTS 

iRobertson et alJ d2006al) andlSoringel & Hernauistl (120051). 
and s ubsequently iGovernato et al.l d2007l) and Naa b et alJ 
(2006) have demonstrated that even major mergers can leave 
remnants with non-negligible disk components. Nevertheless, 
we wish to highlight several properties of these disks first, to 
establish their existence and nature. Moreover, we wish to 
ensure that we can robustly identify disks in our merger rem- 
nants, before going on to analyze the conditions for their sur- 
vival. In order to do this, we have considered several methods, 
include e.g. fitting the surface brightn ess profiles to traditional 
bulge-disk decompos itions, (see e.g. IRobertson et al.l h2006a: 
Hopkins et al. 2008d), kinematic decomp ositions based on 
one and two-dimens ional velocity maps dCox et al.l [2006b; 
iHoffman et alJ[2007h . and three dimensional component fit- 
ting. These ultimately give similar results, although e.g. 
surface brightness profile fits and velocity decompositions 
can be considerably dependent on the viewing angle and 
are not especially robust at separating a small disk (in a 
bulge-dominated system) from e.g. other kinematic subcom- 
ponents or rotat ing bulges (a well-known observational dif- 
ficulty, see e.g. iBalcelis et alJ 120071: lKormendvet"aT] 120081: 
iMarinova & Jo geel 12007; Bara zza et al.ll2.Q08l and references 
therein). 

We therefore choose to take advantage of our full three- 
dimensional information in the simulations to easily decom- 
pose bulges and disks in a simple, automated fashion. For 
convenience, let us consider the remnant in cylindrical coor- 
dinates x = {R, (p) where the axis of symmetry (£) is defined 
by the net angular momentum vector of the baryonic mass in 
the relaxed remnant. The effective rotational support of any 
given stellar or gas particle in the simulation is then 



where v c is the circular velocity 

lGM enc (r) 
vc=]]^— (2) 

(here r is the three dimensional radius from the galaxy center). 
If we consider the distribution of baryonic mass in v ro t, we find 
a clear segregation between bulge and disk components. 

Figure Q] shows this for three simulations with large disks 
in the remnant. There is clearly a bimodal distribution in 
v rot , with one component having relatively little rotation (the 
bulge, with a peak near v lot sa 0), and one component be- 
ing largely rotationally supported (the disk, with a peak near 
^rot ~ 1)- There are two stellar populations in these remnants, 
with a clean division in their rotational support. 

We can, from this plot alone, estimate a robust disk-bulge 
mass ratio, from fitting e.g. the sum of two Gaussian com- 
ponents (disk and bulge) to this distribution (or in a non- 
parametric sense, by assuming the bulge component has a 
symmetric rotation distribution about its peak, and mirroring 
the distribution about that, taking what remains to be disk). 
Our results are not sensitive to the exact details of our de- 
composition, but we experiment with a few different methods 
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FIG. 1 . — Examples of merger remnants with large disks. Top: Edge-on projected stellar surface brightness of the galaxy. Middle: the distribution of all stars 
in their rotational support, v mi /v c (r), where v c (r) is the circular velocity at r and Vrot = v <j> m cylindrical coordinates is the rotational velocity about the net stellar 
angular momentum axis. We decompose the clearly bimodal distributions into bulge (orange, peak near v 10l /v c (r) ~ 0) and disk (blue, peak near v 10l /v c (r) ~ 1) 
components, with stellar mass fraction in the disk component (/disk) labeled. Bottom: Azimufhally averaged face-on surface brightness profile. We show the total 
profile (black) and profile of each of the components separated by their rotational support (orange and blue respectively). We fit the total profile to a standard 
bulge+disk decomposition, and show the resulting fitted bulge (red dotted) and disk (blue dotted) components and disk mass fraction. The two methods recover 
similar decompositions in almost all cases: the "disks" are rotationally supported with extended exponential profiles, the "bulges" are dispersion supported with 
compact Sersic-law profiles. We show three example remnants typical of our simulations: Left: Equal mass (mass ratio 1 : 1) merger remnant with ~ 40% gas at 
the time of merger. The remnant is a bulge dominated elliptical/lenticular, but has a prominent smooth stellar disk with ~ 30—40% of the mass. Center: Major 
(mass ratio 1 : 2) merger remnant with ~ 20% gas at the time of merger. The remnant is a marginally disk-dominated SOa-type galaxy, with some spiral structure 
in the disk. Right: Minor (mass ratio 1 : 8) merger remnant with ~ 15% gas at the time of merger. The remnant is a Sb/Sc disk with a flattened, compact bulge. 



in order to estimate uncertainties on the bulge-disk decom- 
position which we refer to below. Again, we have repeated 
our entire analysis using alternative estimators of the disk-to- 
bulge mass (direct profile fits and velocity profile decomposi- 
tions), and find that the same scalings apply in all cases (the 
uncertainties in the decomposition of a given simulation do, 
however, increase). 

Figure [2] shows one simulation from Figure Q] (SO major- 
merger remnant), using this method to decompose the rem- 
nant into a stellar bulge and stellar disk. We also show the 
gas separately, which can cool and therefore forms an ex- 
tremely thin disk. The properties are exactly what would 
be expected for a typical bulge-disk system: the "bulge" 
is a somewhat flattened ellipse with ellipticity e w 0.3-0.4 
(H/R w 0.6-0.7), and is a pressure-supported system, with 
one-dimensional velocity dispersion a ~ 120- 150kms _1 (de- 
pending on the sightline and slit width) and a rotation veloc- 
ity ~ 30 -50 km s" 1 . The resulting rotation parameter of the 
bulge itself ({V/cr)* ~ 0.4) is typical of reasonably rapidly ro- 
tating bulges. It is also compact (as expected), with projected 
R e ~ 1— 2 kpc. The stellar disk is like that of a combined thin- 
thick disk system, with H/R ~ 0.15-0.2, and exhibits a flat 
rotation curve with V max ^200kms _1 . The disk is rotationally 
supported with typical values for a disk of similar mass and 
overall morphology, V/cr ~ 2-3, and it is far more extended 



than the bulge (R e ~ lOkpc, putting it on the observed disk 
size-mass relation shown below). The properties of the gas 
disk are similar, with the obvious exception that, since the gas 
can cool, it forms a very thin disk (H/R < 0.05). 

Figure [3] shows the components of another galaxy (Sbc 
minor-merger remnant) from Figure [T] in the manner of Fig- 
ure |2 As expected given the smaller bulge-to-disk ratio in 
this case, the system is more flattened, with even larger rota- 
tional support (V/cr ~ 5— 10 and H/R < 0.1 even in the stel- 
lar disk). The bulge is clearly a distinct dispersion-supported 
component, despite being relatively flattened. 

Figure [4] shows the components of the third galaxy (ellip- 
tical major-merger remnant) from Figure Q] in the manner of 
Figure |2] We show this case to demonstrate that embedded 
disks can be recovered reliably, and are indeed real, rotation 
supported (V/cr ~ 3— 5), and relatively thin (H/R < 0.2 in the 
stellar disk, < 0.1 in the gaseous disk) kinematic objects even 
in simulations where they are not a majority of the mass (here, 
we find B/r- 0.7). 

As a further check that these are indeed real disks, Fig- 
ure [5] plots the disks in disk-dominated simulation rem- 
nants on the baryonic Tully-Fisher and stellar size-mass 
rela tions observed for disks of similar morphology (see 
e.g. iBell & de Jondl200l"t lMcGaugrJl2005t IShen eTalTl2003h 
Courteau et al. 2007). For convenience, we will use the esti- 
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FIG. 2. — Bulge (left), stellar disk (middle), and gas (right), in the remnant of a 1 : 2 mass-ratio major merger with ~ 20% gas at the time of the merger, on 
a typical random orbit (center panel in FigurefT}. From top to bottom, panels show 
edge-on scale height H/R of the component as a function of circular radius R = \x\ I 
and velocity dispersion cr(r) (red), (e) Rotational support measure v/cr. 
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FIG. 3. — Bulge, disk, and gas, shown as in Figure|2] for the remnant of a 
1 : 8 mass-ratio minor merger with ~ 15% gas, on a polar orbit. 

mated bulge-to-disk ratio as a proxy for morphology through- 
out, with the values as labeled in Figure [5] We take their 
velocities here from the projected disk rotation curves where 
they are fiat, and take R e as the projected half-mass radius 
(note that this is different from the exponential disk scale 
length h; for a pure exponential disk R e = 1 .678 /z, and we 
convert the observations where necessary accordingly). 

For each morphological class, our simulations agree well 
with the observed Tully-Fisher and size-mass relations. Given 
our limited sampling of very minor mergers with mass ratios 
> 1 : 8, we have only a few simulations with final B/T < 0.2, 
but those nevertheless agree (as expected, since they have only 
been slightly modified from the original disk). We stress that 
we are not claiming to reproduce the Tully-Fisher or stellar 




FIG. 4. — Bulge, disk, and gas, shown as in Figure [5] for the remnant of a 
1 : 1 mass-ratio major merger with ~ 40% gas, on an inclined polar orbit. 

size-mass relation of disks in an a priori manner: our (pre- 
merger) disks are constructed, by design, to more or less lie on 
the observed correlations. What we are saying is that, given 
progenitor disks that are similar to those observed, disks that 
form after or survive mergers (even major mergers) will re- 
main on the appropriate correlations for their stellar mass and 
morphology. In short, when disks do survive mergers, they 
are "real" disks in the observable sense, not highly flattened 
bulges or unusual kinematic subcomponents. 

4. DISK FORMATION IN MAJOR MERGERS 

Clearly, even major mergers can and do produce remnants 
with significant disks. We therefore ask how these disks form, 
and whether we can derive some analytic expectation for their 
masses as a function of progenitor and merger properties. 
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FIG. 5. — Simulated disk-dominated merger remnants on the observed 
baryonic Tully-Fisher relation (top) and size-stellar mass relation (bottom). 
We take V lot from the rotation curves as in Figure [5] and R e is the median 
projected half-mass radius. We compare with th e observed relations as a 
function of morphology from Courteau et al. 1 2007) (solid lines in each panel; 
dotted lines show the observed ±1ct scatter), for SO-Sa (left; red), Sab-Sbc 
(center, green), and Sc-Sd (right; blue and black, respectively) galaxies. For 
convenience, we assign our simulations a "morphology" based on the bulge- 
to-disk ratio as labeled. Most of our simulations are major (mass ratio 1 : 
1) mergers that yield significant, but not dominant disks (and are therefore 
not shown here). Very late types are only produced in our limited subset 
of small mass ratio (> 1 : 8) mergers. In any case, the remnant disks all 
lie on the observed Tully-Fisher and size-mass relations appropriate for their 
morphology - coupled with their rotation and scale heights, we can say they 
are real disks in the observable sense. 

4.1. Components of the Remnant: Surviving Gas Disks 

For simplicity, let us begin with the case of an identical 
1 : 1 mass ratio merger (we will generalize to arbitrary mass 
ratios in § 14.3.31 below). Early in the merger, the galaxies 
experience a first passage and begin to lose angular momen- 
tum to the halos, rapidly coalescing on a timescale of order 
a couple orbital periods. In the final merger and coalescence 
of the galaxies, the stars which are initially "c old" (i.e. pre- 
merg er disks) w ill scatter and v i olentl y relax dLynden-Belll 
1967), forming a ide Vaucouleursl (fl948)-like quasi-spherical, 
dispersion supported profile. Of the gas available at the time 
of the final merger, some will lose its angular momentum, fall 
into the galaxy center, and (given the sudden rapid increase 
in density) rapidly transform into stars in a central starburst, 
forming a compact, central dissipational component of the 
remnant bulge (for a detailed study of this component, see 
iHopkins et alj|2008alldh . Gas that is at sufficiently large radii 
that it cannot efficiently fall in, or gas that for whatever rea- 
son cannot efficiently dissipate or lose angular momentum, 
will rapidly see the central potential relax (the equilibration 
timescale of the central bulge is only ~ 10 s yr) and, having 
conserved its angular momentum, will rapidly cool and re- 
form a thin, rotationally supported disk. iBarnes & Hernquisa 
( 1996) outline this process and show, in detail, how the cool- 
ing gas that survives the merger rapidly settles into a typical, 
rotationally supported exponential disk. This will then form 
stars, which constitute a new stellar disk. 

We emphasize these three components: 

Pre-Merger Stars: These (along with the dark matter) 
constitute the collisionless (dissipationless) component of the 
merger. Because they are collisionless, the stars and dark 
matter distributions mix in the merger. A given star, as it 
moves through the merging galaxies on a random orbit, feels 
a rapidly fluctuating potential, which deflects its orbit and al- 
lows for the phase space distribution of the particles to uni- 
formly mix. This violent relaxation process gives rise to 
a pr essure supported sy stem dominated by random veloci- 
ties dLvnden-Bellll 1 9671) and transforms initially exponential 
disk into quasi-spherical Ide Vaucouleursl (fl948)-like Sersic- 



law profiles. In the limit of a 1 : 1 mass ratio merger, it is 
a good approximation to assume that all of the stars are vi- 
olently relaxed - the merger is sufficiently "violent" that no 
significant component of the pre-merger stellar disks will sur- 
vive the merger. This is not necessarily true at lower mass ra- 
tios (see § 14.3.31 ). but it simplifies our analysis to begin, while 
we consider such mergers. 

The remaining two components of the remnant can be iden- 
tified with the gas supply available at the time of the final 
merger. 

Starburst Stars: This is the remnant of a dissipational star- 
burst, triggered in the merger. Some fraction of the gas will 
efficiently lose its angular momentum in the merger. Because 
gas can dissipate energy, it will then necessarily rapidly fall 
into the center of the merging system (essentially free-falling 
to t he center until the colla psing gas becomes self-gravitating; 
see IHopkins et alj|20 08d). Collecting a large gas supply in 
the center, the result is a rapid, highly concentrated starburst 
- in gas rich cases, this i s analogous to that o bserved in e.g. 
nearby merging ULIRGs jSoifer et al.ll 1984atfbl : IScoville et all 
1986; Sargent et alj I1987L 1 1 989b a n d recent merge r rem- 
nants dKormendv & Sanders! 119921; iHibbard & Yunl 119991; 
Rothberg & Joseph 2004). This builds up a dense, com- 
pact central stellar distribution, that raises the central phase 
space density and yields an effectively smaller, more baryon- 
dominated remnant. The starburst stars, being so concentrated 
(and typically having a more mixed orbital distribution owing 
to the random velocities of infalling gas in the starburst), are 
clearly part of the bulge (although they may have slightly dif- 
ferent Sersic profiles and kinematics from the more extended 
bulge formed from violent relaxation of the pre-merger stars). 
This component is important for the structure and scalings of 
the bulge/spheroid co mponent, and we study it in detail in 
Hopkins et al. (2008a g.d). It is essentially the dissipational 
component of the merger. For our purposes here, however, 
this is the gas "lost," which becomes part of the bulge and no 
longer contributes to the remnant disk. 

Surviving Gas/"Post-Merger" Stars: The gas that does 
not lose its angular momentum will, as described above, form 
a new disk as the remnant relaxes. For a 1 : 1 merger, since (as 
noted above) the entire stellar distribution is violently relaxed, 
the post-merger disks can be entirely identified with gas that 
survives the merger. It is not, in this case, so much that the 
initial disks survive the merger intact, as it is that some of the 
gas remains at large radii/with significant angular momentum, 
which can rapidly re-form the disk after the merger. Essen- 
tially then (for major mergers), the question of how much of a 
disk will remain post-merger is a question of how much of the 
gas (at the time of the merger) will or will not lose its angular 
momentum. 



4.2. How Does the Gas Lose Its Angular Momentum? 

How, then, does the gas lose angular momentum in a 
merger? The basic process has been understood since early 
simulations involving highly simp lif ied models for ga s dis- 
sipation in iNoguchf (1 19871 119881) . iHernquistl d 1989b . and 
IBarne s & Hernquist ( 199 ip7 W ith improved numerical mod- 
els, |Barne!jOlernguist (1996) followed this process in de- 
tail, and showed that what happens in a typical major merger 
is as follows: the non-axisymmetric perturbation (owing to 
the companion) in the system induces (largely after first pas- 
sage and on the final coalescence, since this is where the in- 
teraction is significant) a non-axisymmetric response in the 
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disk. 3 A stellar bar and gas bar form, but because the gas is 
collisional and the stars are collisionless, the stellar bar will 
trail or lag behind the gas bar by a small offset (typically ~a 
few degrees). The stellar bar therefore torques the gas bar, 
draining its angular momentum, and causing the gas to col- 
lapse to the center. 

Figure [6] illustrates this in a couple of representative 1:3 
mass ratio major merger simulations. For a more detailed 
de scription and illustra t ion of the relevant physics, we refer 
to iBarnes & Hernquistl (Qj996) (particularly their Figures 3- 
8); but we briefly outline the scenario here. We show the 
morphology of the gas before the merger, when the disk is 
undisturbed, and shortly after both first and second passages 
(the second passage leading, in these cases, to a rapid coales- 
cence), as well as in the relaxed remnant. 

The bar-like non-axisymmetric perturbation induced by the 
close passages is clearly evident; the stars show a similar mor- 
phology at each time, with a small phase offset in the bar pat- 
te rn and (in the remnant) a stellar bulge. As we discuss in 
§ 14.3.21 a prograde encounter, being in resonance, produces 
a noticeably more pronounced bar distortion (both in ampli- 
tude - effectively "bar mass" - and spatial extent). Shortly 
after each passage, this double bar system efficiently removes 
angular momentum from the gas, allowing it to fall into the 
center of the galaxy and participate in a centrally concentrated 
starburst. 

Following Ba rnes & Hernquistl (Q996), we track the gas in 
the primary disk that will turn into stars in the final starburst, 
calculating the net instantaneous gravitational torque deceler- 
ating the disk rotation. We can coarsely infer what the total 
effective torque must be by simply differentiating the specific 
angular momentum of this gas at a given time, and compare 
this to the net torque from different sources. Specifically, we 
separate the instantaneous gravitational torques into the inter- 
nal torques - those from the stellar disk in the same galaxy as 
the gas, chiefly from the bar (since the axisymmetric disk, by 
definition, exerts no net torque) - and the external torques - 
those from the gravity of the secondary galaxy itself and the 
extended halos and their substructure. 

It is clear that, especially for the phases of interest shortly 
after second passage and leading into the final starburst, when 
this gas loses its angular momentum, the total torques are 
dominated by internal torques from the stellar disk/bar sys- 
tem. The agreement between these torques and the rate 
of change in the specific angular momentum further argues 
that there are no other major sources of angular momentum 
loss (specifically, both this comparison and direct calculation 
demonstrate that the "hydrodynamic torques" defined by pres- 
sure forces are not dominant). 

As a result of these torques, gas within some critical radius 
where the internal torques are strong (roughly inside the "bar 
radius" in Figure[6]) rapidly loses angular momentum. We de- 
fine this radius more precisely in § l4.3l below, but it is clear 
in the figure that at sufficiently large radius the bar perturba- 
tion is weaker (and moreover, at larger radius the potential of 

3 In what follows, we will refer to this non-axisymmetric response as a 
"bar," for simplicity and because morphologically the induced feature, at 
least for some time during the merger resembles bars in isolated barred spi- 
rals. However, we caution that the formation mechanism which excites this 
response may be different from that causing bars in isolated galaxies. Fur- 
thermore, while the non-axisymmetry is present throughout the merger, it at 
times would not be classified as a bar morphologically, particularly during 
the final coalescence of the galaxy nuclei, when the resulting gas inflows are 
strong. 



the disk, whether barred or unbarred, appears increasingly ax- 
isymmetric); gas outside of these radii is relatively unaffected. 
In general, then, the means for a more efficient encounter to 
consume a larger fraction of the gas in the disk is to induce a 
stronger bar disturbance, which is able to effectively exert in- 
ternal torques out to larger radii, stripping more gas of angular 
momentum and bringing it into the central starburst (as evi- 
dent in the stronger prograde encounter in Figure [6). Finally, 
the system relaxes - the gas that has not been subjected to 
strong internal torques, having retained its angular momentum 
(at least in large part), can rapidly re-form a disk. This may 
entail some redistribution of that angular momentum ("filling 
in" where the bar depleted the gas of the disk), but does not 
le ad to further significant angu lar m omentum loss. 

IBarnes & Hernquistl d 1996b and IBarnes! d!998h illustrate 
that this internal torquing is by far the dominant source of an- 
gular momentum loss, for typical orbits. This is because the 
stellar bar is: (a) more or less aligned in the plane with the gas 
bar, (b) trailing it by a small amount, and (c) relatively long- 
lived (it lives the rest of the duration of the merger, as opposed 
to the short time that is e.g. pericentric passage). The compan- 
ion itself (either its baryonic mass or its halo), in most orbits, 
is not perfectly aligned with the gas disk, and the torque di- 
rectly from it is much weaker (the tidal torquing drops by a 
factor ~ (^disk/^peri) 3 )), and it can act only for a short du- 
ration on pericentric passage. There are some pathological 
orbits (e.g. perfectly coplanar prograde orbits) where this is 
not true, but these are exceptional cases, and we discuss them 
in§EL6] 

At the final merger, one might image that mixing of random 
gas orbits or collisions and shocks would rapidly drain angu- 
lar momentum, similar to what happens to the stars in violent 
relaxation. However, this is not possible, precisely because 
the gas is collisional: a Lagrangian gas element cannot go 
back and forth through the galaxy, but sticks to the other gas 
which has some net angular momentum. There could in prin- 
ciple be some net angular momentum cancellation, but this is 
inefficient - the net angular momentum will almost always be 
comparable to the initial total. Even assuming random cancel- 
lation between two disks with comparable absolute angular 
momentum, the average change in net specific angular mo- 
mentum is a factor ^2/3; when one accounts for the angular 
momentum of the orbit - typically comparable or even larger 
than that in the disks - there is often no change or even a net 
gain in the gas specific angular momentum in a merger. 

A proper calculation shows that over the range in mass 
ratios fi ~ 0.1 — 1, for a range of typical impact parameters 
b ~ 0.5- 5 Rd, the expected final specific angular momentum 
after cancellation is approximately equal to the initial specific 
angular momentum of the primary (with ~ 20% scatter). Can- 
cellation is therefore inefficient. Even these cancellations, we 
find in detail, do not generally lead to a starburst in the same 
manner as a merger-induced bar, but simply lead to moderate 
disk contraction (and an equal number of mergers will scatter 
towards the opposite sense leading to disk expansion, keeping 
a mean specific angular momentum that is constant). They 
do not cause a starburst because, if two random parcels or 
streams of gas shock and lose angular momentum, the align- 
ment and relative momenta would have to be near-perfect for 
them to lose, say 95% of their angular momentum and fall 
all the way to the central ~ lOOpc where a nuclear starburst 
would occur. Rather, they will lose some fraction of order 
unity of their angular momentum, fall in to a slightly smaller 
radius, and continue to orbit. 
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FIG. 6. — Illustration of the key processes that drive starbursts in a merger. Top: Projected gas density (as in Figureff} in the plane of the disk at representative 
times in a retrograde 1:3 merger (left to right: before interaction, just after first passage, just after second passage/coalescence, after relaxation). For clarity, just 
the gas from the primary is shown. Middle: Same, for a prograde encounter. The passage of the secondary induces a bar-like non-axisymmetric disturbance in 
the primary, which survives after the short-lived passage and removes angular momentum from the gas, leading to a starburst. The same process occurs on both 
passages, with a larger (albeit less "bar-like") asymmetry on coalescence. The prograde encounters, being in resonance, induce a stronger response that extends to 
larger radii. Bottom: Quantities of interest in the merger. Left: Star formation rate as a function of time in the prograde encounter (solid line; retrograde is similar, 
but with a weaker enhancement in bursts). Red dotted line shows the specific angular momentum of the gas that will participate in either burst, in arbitrary units); 
as the gas rapidly loses angular momentum after the passages, it drives a central starburst. Center Left: Net specific torque on the primary gas that will participate 
in the final, central starburst, as a function of time in the merger (in units of the initial total angular momentum per Gyr). We compare the roughly numerically 
estimated net torque (diamonds; from differentiating the specific angular momentum of the gas) and the torque from two sources: stars in the same disk as the 
gas (internal torques; black thick line), and the secondary galaxy and extended halos (external torques; red dot-dashed line). The loss of angular momentum that 
drives the secondary burst at t ~ 2.2-2.5 is driven by internal torques from the disturbed stellar disk; not the torque from the secondary galaxy itself. Center 
Right: Final specific angular momentum content of material (gas plus stars) that was originally gas at R < R alt (our predicted radius where merger-induced 
internal torques should be efficient at removing angular momentum) or at R > R cr \t. There is a strong division: gas inside a characteristic radius (corresponding 
to where the internal asymmetry is strong; akin to the co-rotation resonance) is mostly stripped of angular momentum. Gas at larger radii conserves sufficient 
angular momentum to maintain a disk at similar R e and V c . Right: Original (cumulative) radial distribution of gas that participates in the final starburst, relative 
to the initial disk effective radius (same in all cases) for the prograde and retrograde cases shown and a more minor retrograde merger. More resonant (prograde) 
and more major encounters induce a stronger response, with a larger co-rotation radius, and so torque gas out to larger radii and efficiently strip more gas of 
angular momentum as predicted. 



Without the bar that can continuously drain angular mo- 
mentum, the true burst is indeed inefficient. This ini- 
tial bar-induced angular momentum loss scenario has been 
well-establishe d in subsequent numerical studies (see e.g. 
Noguchil 119871 119881: iHernauistj [l989t iHernquist & Barnesl 
1 99 U iBorderies et al.ll 19891: iBarnes & Hernquistil 1 99 U 1 1 9961: 
Barnes 1998; Mihos & He rnquistil 996tlSpringel & Hernauistl 
2005b iRobertson et al I l2006at ICoxetalJ l2006bl l2008at 
Berentzen et al.l 2007; N aab et al.ll2006l) . We therefore can 
simplify our question to ask: how efficient will a given lag- 
ging stellar bar be at removing angular momentum from a 
leading gas bar? 

4.3. A Simple Model: Overview 

Consider a disk that contains a total gravitational mass M 
(which can include a bulge and dark matter as well; the disk 
mass fraction we will denote /disk) within a characteristic 
scale length Some convenient dimensional variables are: 

M disk = /diskM (disk mass) 



M Dar = fbaiM (stellar bar mass) 



GM 



v c = \ (characteristic circular velocity) 

V ^d 

Vc 

fld= — (characteristic frequency) 

^d 

P = — (rotation period) . 



(3) 



We also define the disk thickness according to a characteris- 
tic (assumed constant) relative scale height (height H versus 
radius R; H/R =constant) 



H = H/R (disk scale height). 



(4) 



In these units, the circular velocity at a given cylindrical ra- 
dius R is given by 



VcircM = v c v(r) = v c 



M enc (r) R d 
R 



M 



(5) 
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and dimensionless lengths are defined by 

x = x/Ra 
y = y/R d 



R = Vx*+f/R d 



(6) 



Throughout, we will use this notation: e.g. the dimensional 
variable u is equal to the dimensionless variable u times the 
appropriate combination of dimensional constants above. 

We will show that, in such a disk, a gas bar with a lagging 
stellar bar will efficiently cause gas to lose its angular mo- 
mentum and dissipate into a central starburst. This will be the 
case for gas interior to a radius 



< a(l-f gas )f disk F(9,b)G(ji), 

Rd 



(7) 



where a ~ 1 is an appropriate integral constant (depending 
weakly on details of the stellar profile shape and bar dynam- 
ics), /g as is the gas fraction in the disk and /disk is the disk 
mass fraction. The factor 



GO) = 



2/i 
(1+M) 



(8) 



contains the dependence on the merger mass ratio (where \i < 
1 =M 2 /Mi). The term 



F(9,b) 



1 



i + [b/R d y 



3/2 



1 



(9) 



accounts for the orbital parameters: b is distance of pericentric 
passage on the relevant final passage before coalescence (~ 
1-a couple Rd, for typical cosmological mergers) and f2 is 
the orbital frequency at pericentric passage, 



Vperi Rd ,~ 

— cos(#) 



: v/2(l+Ai)[l + (^d) 2 ]" 3/4 cos«?) 
;0.6 cos(6»), 



(10) 



where 8 is the inclination of the orbit relative to the disk, and 
the last equality comes from adopting typical cosmological 
orbits and mass ratios (but in any case, this is quite weakly 
dependent on the mass ratio). 

In the following sections, we derive this scaling piece by 
piece, and compare each aspect to the results from our library 
of hydrodynamic simulations. We show that it is robust and 
accurate as an approximation to the behavior in full numeri- 
cal hydrodynamic experiments over a wide dynamic range of 
several orders of magnitude in surviving disk fraction (from 
systems with ~ 80- 100% of their disks surviving a merger to 
systems with < 1 % disk after a merger), as well as the entire 
dynamic range in mass, gas content, orbital properties, and 
different feedback prescriptions with which we experiment. 

4.3.1. Dependence on Disk Gas Content 

Let us consider an infinitely thin gas bar (a good approxi- 
mation, owing to the efficiency of gas cooling) in a potential 
that is otherwise cylindrically symmetric except for the pres- 
ence of a stellar bar. For simplicity, assume that the gas bar 
follows a fixed pattern speed fib (~ fi ; we will derive the 
pattern speed later) in the disk (while there is not exactly a 
constant pattern speed in the outer regions of the disk, the 
torques are weak there in either case, and this approximation 
is globally quite good). We take z = to be the plane of the 



disk, and, without loss of generality, consider a frame rotating 
with the pattern speed of the gas bar, so that the bar lies along 
the x axis. The material in the bar is rotationally supported, 
so it has instantaneous velocity = — v c v(R)y, where v c v(R) 
(defined above) is the circular velocity at each point x. 

Now, consider a stellar bar of total mass Mbar also at fixed 
pattern speed, but offset by some instantaneous angle fa from 
the gas bar (i.e. along the axis y = tan(0/,)x). The mass 
per unit length in the bar at a distance R along the bar is 
dM bar /d/? = (Mb ar /^d) ^(R/R d ), where £ is the appropriate di- 
mensionless mass profile and Rd is some characteristic scale 
length (usually corresponding closely to the scale length of 
the unperturbed disk). If the initial disk is in equilibrium, (i.e. 
if the bar is some reasonable perturbation to the initial sys- 
tem), then the unperturbed net acceleration in the x direction 
at some point x in the gas bar will just be cancelled by the 
rotation of the system. Of interest here is the torque; if the 
stellar bar is also thin, then at a point x = xRd in the gas bar, 
the net torque per unit mass from the stellar bar will be 



d J ~ D dv y 
— = xRd — - 
df dt 



GM ba 
Rd 



k(fa,x), 



(11) 



where Iq ~ 1 is a dimensionless integral which depends 
weakly on fa and x (at large fa, Iq — > 0, reflecting the fact 
that the torque is dominated by times when the bars are close; 
since <§; 1 is expected, it is a good approximation to ignore 
the fa dependence of 7o). 

If we assume the stellar bar is infinitely thin, there is a weak 
divergence in Iq as fa — > (the accelerations become large 
when the bars nearly overlap). More accurately, we can allow 
for some finite height in the stellar disk/bar (it will always be 
thicker than the gas disk/bar). Let the stellar bar have a con- 
stant relative scale height H/R given by Equation ©, and for 
simplicity take its vertical profile to be constant density out 
to a height ±H (although assuming a more realistic vertical 
profile oc exp(—\z\/H) or oc sech 2 (z/i/) makes almost no dif- 
ference to our calculation). The specific torque at x in the gas 
bar now becomes 



dt 



GM b!l 



Rd y/ 



1 



sin" 



h(fa,x,H), 



(12) 



where 7) is an even weaker function of fa and x than Iq. The 
important behavior is entirely captured ignoring h, namely 
that the finite width of the stellar bar suppresses the numerical 
divergence seen earlier. 

The bar mass Mb ar represents the stellar mass in the disk 
that is effectively part of the bar at the appropriate instant. We 
can therefore parameterize / bal - = M bal -/M as 



/bar = (1 "/gas) /disk ^b 



bar ' 



(13) 



Here, /d; s k is the disk mass fraction, and /g as is the gas frac- 
tion in the disk (since we are interested in the cold gas, we 
explicitly ignore gas in e.g. a bulge or hot halo component). 
Therefore, the stellar mass of the disk is (1 -/ gas )/disk - this 
defines the maximum mass that could be in the stellar bar. 
The parameter ty' bgx thus defines the bar "efficiency" - in an 
instantaneous sense as we have defined it here, ^ bal . = means 
there is no stellar bar, ty' bai = 1 implies the maximal stellar bar. 

Already, we have one significant scaling - the bar strength, 
and correspondingly the strength of the torques on the gas, 
scale with (1 -/ gas )- In very gas rich systems where / gas — * 1, 
there is no stellar mass to form a lagging bar and remove an- 
gular momentum from the gas. The gas itself may form a bar, 
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but without a stellar bar to drag it, the angular momentum 
loss (over the timescales of relevance for a merger 4 ) is inef- 
ficient. This is well kno wn in e.g. dynamical studies of pure 
gas and stellar bars (e . g. | Schwarzl [19811: lAth anassou la et alJ 
19831: iMennigen ll984t1Combes et all Il990t iFriedli & Bend 
1993b lO'Neill & Dubinskill2003l) . There might be some an- 
gular mo mentum loss in such a case, b etween e.g. bar and 
halo (e.g. lHernquist & Weinberd [l992a). but it will be small 
- certainly nowhere near the efficient stripping of angular mo- 
mentum needed to induce a significant starburst. 

If we consider a Lagrangian gas element at some initial ra- 
dius xo, then its orbit will decay as it loses angular momen- 
tum. The instantaneous rate of change in the radius R = \x\ 
will be given by dR/dt = (R/v^dv^/dt. The characteristic 
timescale for the system to evolve is given by 27r/f2d, where 
~ v c/h is the characteristic frequency of the disk. Define 
the timescale 



n d v c GM 
2tt R d \ Rl 



(14) 



where M is the total effective gravitational mass of the disk 
and R d is again a characteristic scale length. We now have: 

di 2tt 1 GM bar 

^-=o -5 — 1\ = 2lT fbaxhwb-, X) . (15) 

dr S2 d 

Because the merger occurs on a couple of dynamical 
timescales, i.e. a time At ~ 1, to lowest order (ignoring e.g. 
the complications of different orbital parameters) we expect 
that gas within a radius 



x < At 



di 



(16) 



will efficiently lose angular momentum and fall to the center, 
becoming part of the central starburst, while gas at I ^ At g| 
will avoid the starburst. This defines a scale 



^1 <n. 
Rd ~ 



/ga S )/disk*bar(^,H,...) 



(17) 



within which the gas will lose angular momentum. For con- 
venience, we have collected all of the dimensionless integral 
factors, including $ bal (the efficiency of forming the stel- 
lar bar) and the dynamical integral factors (e.g. h, h) from 
above, into the term ^bar that represents the full solution. We 
write v I'bar(0fc, H...) because, as we will show, this quantity (at 
present) encapsulates our ignorance of e.g. the orbital param- 
eters and merger mass ratio; for a 1 : 1 merger on a typical 
orbit, however, ^bar ~ 1 ■ 

The total gas mass which will lose angular momentum and 
fall into the center of the galaxy will be / gas x /(< 7? gas ), where 
/(< R gas ) is the mass fraction within the characteristic radius 
above, according to the details of the mass profiles and dimen- 
sionless integral above. We consider solutions for a variety of 
p rofiles, includin g e.g. an exponential, isothermal sphere, and 
a iMestell (1 19631) 1 /R disk profile. In general, we find that 
there is little difference between the predictions for these var- 
ious profiles - the differences in the mass profile shapes tend 
to cancel out and leave only weak corrections to the simple di- 
mensional scaling. An exponential disk with £ oc exp(-R / R d ) 
contains a mass fraction 



l-(l+/?//?d)exp(-/?//? d ), 



(18) 



At least, in this case, a major merger. The situation becomes m ore com - 
plicated in the limit of minor mergers with mass ratios ~ 1 : 1 0; see § 14.3.71 
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FIG. 7. — Mass fraction formed in the central, dissipational starburst as a 
function of gas mass fraction at the time just before the starburst, in a suite 
of major 1 : 1 mass-ratio mergers. Solid l ines are our theoretical predic- 
tions (/burst = /gas(l — /gas)^, see Equation 1 19k dotted line corresponds to 
bursting all the available gas (/i, urst = /gas). We show results here for several 
orbits from TablefT] a typical random orbit with both disks inclined (e: top), 
an inclined polar-prograde orbit (k: middle), and a polar-polar orbit (f: bot- 
tom). The simulations agree well with our analytic predictions: more gas-rich 
mergers are less efficient at torquing angular momentum away from the gas 
and funneling it into the starburst (efficiency ~ (1 — /gas))- 

within a radius R; here 'J'bar must be solved numerically. We 
obtain nearly identical predictions, however, assuming a l/R 
disk or an isothermal sphere profile for the gas, which allows 
us to analytically solve the relevant equations and write the 
predicted gas fraction consumed in the form: 

/burst = /gas ( 1 - /gas) /disk *bar(<&> ,H,...), (19) 

where 'J'bai ~ 1 can be analytically calculated for these pro- 
files (with the equations above) under certain conditions: if 
the dependence on the orbital parameters is separable and we 
define ^bar by the requirement that the radius in Equation dl7l i 
satisfy x = At then for instantaneous bar lag of (pi, ~a few 
degrees in a thick disk of height H ~ 0.2, ^bar is given by 

sin(2 <f> h ) 



;F(. 



»{ 



exp 



sin 2 (0 fe ) + i/ 2 



1. 



(20) 



where we explicitly show F(...) as we have suppressed our 
ignorance of the orbital parameters. Nevertheless, this simple 
scaling alone provides a remarkably successful description of 
many of our simulations. 

Figure [7] tests this simple prediction. For a suite of merger 
simulations, we compare the mass fraction in the central star- 
burst, /b ul st, to the gas content of the (immediately pre-merger) 
disks, / gas . We can either determine the starburst mass frac- 
tion by directly measuring the gas mass that loses its angular 
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FIG. 8. — Relaxed post-merger remnant disk mass fraction versus gas frac- 
tion just before the merger, for 1 : 1 major mass-ratio mergers. In this case 
essentially all the pre-merger stellar mass is transformed (violently relaxed) 
into bulge - the disk is formed from the gas that survives the merger. Panels 
consider different orbits, with points as Figure|7] Solid line s ar e our theoret- 
ical predictions (/disk = /gas[l — (1 — /gas) see Equation ll9t . dotted lines 
correspond to all the gas surviving and forming a disk (/disk = /gas). Again, 
the simulations agree well with our analytic predictions; gas-rich mergers are 
inefficient at stripping angular momentum from the gas, leaving significant 
gas content that rapidly re-forms a post-merger disk. 

momentum and participates in the brief nuclear starburst, or 
by measuring the gas content that survives and forms a disk 
(described in § and assuming the gas that did not survive 
(relative to that available just before the final merger) was part 
of the burst. In either case, we obtain a nearly identical answer 
for each simulation. For now, we consider only simulations 
with a 1 : 1 mass ratio - we will generalize to arbitrary mass 
ratios below. In all these simulations, the pre-merger /di s k ~ 1, 
and measuring 4>b an d H just before the merger we expect 
^bar ~ 1. We consider one set of orbits at a time - i.e. com- 
pare only systems with the same orbital parameters, so that we 
can temporarily suppress the dependence on them (this yields 
a systematic offset between each set of orbital parameters - 
the solutions plotted account for that following our solution 
in the next section). At a fixed orbit, for these mergers, then, 
the only parameter that matters should be / gas . The simula- 
tions at each orbit span a wide range in / gas , from ~ 0.01 - 1. 

We compare the relation between /burst an d /gas resulting 
from the full numerical experiments to the simple scalings 
predicted by Equations (fT7li-(fT9]l. In detail, we show two 
solutions - first, the scaling given by Equation ( TT9l >, fhi irst oc 
/ gas (l -/ gas ), appropriate for an isothermal sphere or iMestell 
( 1963) disk profile with $bai- ~ 1; and second, the appropriate 
numerical solution (following Equations [T7l[T8T > for an expo- 
nential disk. In either case the analytic solutions are similar, 
and agree well with the trend seen in the simulations. It is 
clear that the efficiency of the burst in simulations is - as we 
predict - not constant. It is not the case that the entire gas sup- 
ply is always stripped of angular momentum and consumed in 
the final merger (which would yield /b ust = / ga s)- Rather, when 
/gas is sufficiently high, only a fraction ~ (1 — / gas ) of the avail- 
able gas is able to efficiently lose its angular momentum and 
participate in the starburst. 

Figure |8] repeats this comparison in terms of the surviving 
disk mass. We argued that the gas that does not lose angular 
momentum in the merger will survive to re-form a disk. Be- 
cause these are 1 : 1 mergers where we can safely assume the 
entire stellar disks are destroyed, we expect then that the disk 
mass fraction will be /disk = /gas - /burst- Using the method de- 
scribed in §|4]to estimate the remnant disk mass fractions, we 




FIG. 9. — Relaxed post-merger remnant disk mass fraction versus our an- 
alytic predictions as a function of gas fraction and orbital parameters, for 
1 : 1 mass-ratio mergers. Error bars correspond to variation using different 
methods to estimate the disk-bulge decomposition. 

plot /di s k versus / gas for each of several orbital parameter sets. 
Again, the exact details of the predictions depend on orbital 
parameters in a manner we derive below, but for now we are 
interested in whether or not they obey the predicted scaling 
with / gas . Indeed, they do. Over 2-3 orders of magnitude 
in fractional disk mass (and ~ 5 - 6 in absolute disk mass), 
the simple scaling here agrees well with full numerical exper- 
iments. It is clear that some of / gas is consumed, as expected 
(if all the gas survived, we would obtain /di s k = / gas ; but in 
fact, especially at low / gas , the efficiency of angular momen- 
tum loss is high as predicted and the gas participates in the 
starburst). 

Figure [9] simplifies this - we again compare /j; s k and / gas , 
but effectively put all the orbits on the same footing by 
plotting /di s k versus the predicted /disk(/gas, •■•) (i-e. includ- 
ing the orbital parameters according to the predictions in 
§ 14.3.21 below). Essentially this amounts to implicitly in- 
cluding F(...) in Equation (fT9l above. The remaining scal- 
ing should just represent the predicted /burst x /gas(l ~/gas)- 
For each simulation, we show an error bar corresponding 
to the range of /disk estimated using different methods (e.g. 
a full three-dimensional kinematic decomposition, one and 
two-dimensional kinematic modeling, and surface brightness 
profile fits, as described in §|4|i. The agreement is surprisingly 
good, given the simplicity of our derivation. Moreover, the 
scatter is quite small - a factor ~ 2-3 at very low / gas and con- 
siderably smaller (< 50%) at high / gas . It seems that our sim- 
ple scaling indeed captures the most important physics of an- 
gular momentum loss - namely that with less fractional stellar 
material in the disk, there is less mass available to torque on 
the gas bar in a merger, therefore less angular momentum loss 
in the gas. 

4.3.2. Dependence on Orbital Parameters 

We now turn to how the details of the orbit affect the loss 
of angular momentum. Before, we made the simplifying as- 
sumptions that the stellar bar lagged by some constant angle 
4>i, and that the characteristic time for the perturbation to act 
was of order the disk rotational period/dynamical time. While 
these turn out to give reasonable scalings, we can improve 
upon them. 

The secondary ("perturber") galaxy will have some charac- 
teristic orbital frequency fi , approximately given by 

n ~^, (2i) 

b 



where v per i is the velocity at pericentric passage and b is the 
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impact parameter or pericentric passage distance. Because the 
behavior we are interested in is relatively short-lived, this is 
reasonable even for first passages or "flyby" encounters - we 
are interested in the orbital frequency at pericentric passage 
because this is when the forcing is strongest and the bar is 
driven. 

Now, consider the frame rotating with the disk/bar at a fre- 
quency ~ fid- In this frame, the secondary galaxy will have an 
apparent frequency for an orbit projected into the disk plane 
of 

fi e ff = Sl cos 9 - fid j (22) 

where 9 is the inclination of the orbit relative to the plane of 
the disk (in standard parlance for orbital parameters as de- 
scribed in §|2]i. Note that we are interested in the component 
of the orbital motion in the plane of the disk: in terms of the 
standard orbital parameters 9 and tj> of the primary galaxy (an- 
gle of the angular momentum vector of the disk relative to the 
plane of the orbit), this is f2 cos 9. For a case with e.g. 9 = 
(prograde) and a parabolic orbit with small impact parameter 
(£1 ~ S7dX the system is maximally prograde - in the frame 
of the rotating disk there is almost no net circular motion of 
the secondary. For the same orbit but 9 = 180° (retrograde), 
the secondary completes a circular orbit around the disk in 
just half a disk dynamical time (f2 e ff ~ -2 fid). The time re- 
quired for the secondary to complete a revolution in this frame 
is therefore (in our dimensionless units r = f/(27r/fid)) 



1 



'cue — q 



(23) 



The timescale for a gas element to lose its angular momen- 
tum and fall to the center of the galaxy is given by our earlier 
estimate of the torque, as ti oss ~ x/(dx/dr). If ti oss <C T c ; rc 
at a given radius, then the derivation we have obtained is es- 
sentially valid: the system sees a quasi-static perturbation to 
the potential, loses its angular momentum, and collapses be- 
fore the perturbation can damp out or circularize. However, 
if ti oss ^> Tcirc, then the system has not lost much angular mo- 
mentum by the time the secondary completes a revolution, 
and will gain some of those losses back as the system comes 
around the other side. In the limit where T a kc is short (much 
shorter than the local dynamical time), for example, then the 
potential is effectively circularized - the gas at these radii may 
undergo oscillatory motion and even have e.g. spiral waves 
driven by this external forcing, but there is no means by which 
the system can introduce a strong net asymmetry to drive in- 
flows. 

We can therefore improve our previous estimate: instead of 
taking At ~ 1 (i.e. a disk rotation period) as the only charac- 
teristic timescale, we argue that gas with 



di/d7 



< 



1 



l-^cos6> 



(24) 



will lose its angular momentum, while gas at larger radii will 
not. 

In detail, we can integrate the equations from I4.3.T1 for a 
parcel of gas at some initial radius xq, in a time-dependent 
potential of this nature. For simplicity, we assume that the 
secondary drives a circular perturbation in the potential with 
frequency uo = Vl Q and calculate the bar response using the 
gaseous disk (assuming, again, that it is infinitely cold) and 
stellar disk (assuming that the scale height H translates into a 
correspondi ng velocity dispersion a/v c ) wave dispersion re- 
lations from Bi nnev & Tremairiel d!987l) . In practice, we find 



this is not much different from assuming that the lag in the 
stellar bar grows with time oc r c i rc (i.e. that the stellar bar can 
keep up or reverse sense tracking the perturbation without sig- 
nificant energy loss; or, more or less equivalently, that the two 
bars are only in phase when the perturbation is strong, and 
then rapidly fall out of phase - at the > 5 - 10° level, once the 
perturbation is weak or reverses its sense as the phase of the 
secondary reaches > 7r/2). In principle, we now have a phys- 
ically motived and fully time-dependent model for 4>t,(t) and 
the response of the gas bar. This allows us to properly inte- 
grate out the dependence of ^bar on and instantaneous con- 
ditions and replace it with the appropriate integral dependence 
on orbital parameters and disk gas content and structure. 

We find that there is a strong division in expected behav- 
ior, at more or less exactly the characteristic radius implied 
by Equation (l24l i. Within this radius, gas (in our simple nu- 
merical calculations) is effectively torqued efficiently as it en- 
ters the gas bar near resonance (but slightly leading the stellar 
bar), and plunges to the center. Gas outside this radius begins 
to feel a perturbation, but then the phase of the secondary cy- 
cles around and the sense of the torques begin to weaken or 
reverse (depending on the details of the orbit), and generate 
wave motion in the gas but no significant angular momentum 
loss or infall. Not only is the transition between these two 
regimes predicted by the simple scalings above, but we find 
in more detailed numerical calculations that the width of the 
transition region (where behavior is more sensitive to the de- 
tails of e.g. the profile shapes and assumptions about the bars) 
is quite narrow, ~ 20% of that radius. 

This should not be surprising. Essentially, what we have 
derived is a rough equivalent of the co-rotation condition, but 
for forced bars as opposed to isolated self-generating (swing- 
amplified) bar instabilities. It is well-known from studies 
of idealized ba r s (see e.g. ISchwarz 1 198 It iPfennigerl 119841: 
i Noguchil 119881; iBinnev & TremaindTT987t iBerentzen et alj 
120071) that gas can be efficiently torqued inwards inside of 
the co-rotation resonance (in the language above, given the 
forcing with pattern speed f2 , this is interior to the radius 
where the relative motion of the secondary is slow relative to 
the dynamical time, and so the perturbation does not circu- 
larize). Moreover, the resonant structure around these radii 
is known to be shar p; if we follow a derivation similar to 
iBorderies et al.l (Il989h (their derivation is intended to apply to 
planetary disks with satellites, but the relevant physics is sim- 
ilar) it is straightforward to show that the detailed numerical 
prefactors will be swamped for all but a narrow range of radii 
around this resonance by the strong dependence of the reso- 
nant forcing on radius (roughly going as some large power of 
(r/r^it) - such that the forcing is strong inside the resonance 
and rapidly weakens outside). 

This gives us confidence that we can adopt the scalings 
above and robustly assume that there is indeed a character- 
istic radius (depending in detail upon orbital parameters) in- 
terior to which the gas will lose its angular momentum. This 
resonant structure of the angular momentum loss is actually 
quite convenient from an analytical perspective, as it means 
that more subtle issues of e.g. the thermal pressure and state 
of the ISM, stellar and AGN feedback, and the exact mix of 
e.g. gas and stars or density structure of the gas will not con- 
tribute significantly to determining which gas can or cannot 
lose angular momentum. Unlike e.g. a self-generating bar in 
an unstable disk, there is no issue of stability analysis - the 
torques inside this critical radius (and the inducing perturba- 
tion) are sufficiently strong such that all the material therein 



How Do Disks Survive Mergers? 



15 



loses angular momentum in a very short time (much less than 
a single orbital time, in practice). 

For example, it is well known that in isolated cases, a 
pure gas disk is more unsta ble to gravitational perturb a- 
tion than a stellar disk (s ee e.g.lChristodoul ou et al.lfl995alfbT : 
Mayer & Wadsley 2004), however in the driven case this is 
not applicable: the distortion in the local stellar/gas distribu- 
tion is caused by the secondary, not by e.g. orbital "pileup" 
or instability in the primary. The location of the resonant ra- 
dius is not determined by the internal structure of the primary 
(unlike in an isolated case, where it is determined by how e.g. 
those orbits can overlap and where various stability criteria 
are satisfied), but rather by the orbital motion of the secondary 
(relative to the internal motion of the primary), and therefore 
knows nothing about e.g. the gas to stars ratio, phase structure, 
and feedback situation in the primary. Inside this radius, the 
distortion is sufficiently strong that it does not matter whether 
one configuration or another is more or less prone to gravi- 
tational instability - the driving force (and therefore angular 
momentum loss) is large in any case. 

Exactly what the pressure support of the gas inside this ra- 
dius is may effect e.g. how far it free-falls after losing angular 
momentum before shocking and forming a central starburst, 
but it will not change the fact that the angular momentum loss 
is efficient. Quantitatively, the torque is 3> j'disk &d (as it must 
be in order for the gas to lose its angular momentum in much 
less than an orbital period); but e.g. the pressure gradients re- 
sisting gas collapse cannot be larger than (in energetic terms) 
HMdV] ~ H 7di S kf2d <C jdisk^rf (or else the disk could not 
be thin) - therefore whether or not there is even considerable 
pressure support or e.g. thermal feedback or a modified ISM 
equation of state makes a negligible correction to the behavior 
seen in the simulations. 

Before moving on, we would like to translate the gen- 
eral scaling above in terms of £1 and Slj into more conve- 
nient parameters. As noted above, il Q ~ v per i/b. We expect 
b ~ 1 - 3 Rd for common parabolic cosmological orbits - as 
we discuss below, for orbits with larger b that will eventu- 
ally merge, all that matters in terms of the end product is the 
impact parameter of the final passage or two when the most 
dramatic forcing occurs, so even for initially larger passages, 
angular momentum transfer to the halo will ensure a value in 
this range towards the final stages of the merger. 

Assuming a parabolic orbit, v pe ri will be given by the in- 
fall velocity from infinity, ^/(GM{1 + [i)/b) (where fi is the 
merger mass ratio, discussed below). Because the merging 
systems are extended, as b — > these expressions should be 
replaced by a more complicated function of b/R& (for the case 
b = 0, the infall velocity asymptotes to the escape velocity 
from the center of the primary ~ -J GM/R&), which requires 
a numerical solution for an arbitrary density profile. In prac- 
tice we find that we can interpolate between the limits b = 

and b Ra quite accurately by replacing b with yj b 2 + R$ 

(which also happens to be an exact solution for e.g. a Plum- 
mer sphere density profile). Combining these factors, we find 
that (for the regime of typical interest) 

7~ COS (6) 

S2 d v c b 

= ^2{\ + ^)[\ + (b/R A fr V4 cos(0) 
« 0.6 cos ((9), (25) 
where the last term comes from inserting a typical major 



merger mass ratio and b ~ 2Rd- The orbital dependence is 
then - as we would expect - largely a function of the inclina- 
tion angle 8. Prograde orbits induce a strong bar response 
- despite the fact that in these mergers the orbital angular 
momenta are all aligned, we actually expect the most angu- 
lar momentum loss and least efficient disk formation. Ret- 
rograde and polar mergers, on the other hand, despite having 
completely un-aligned or cancelling total angular momentum, 
should most efficiently form disks. 

Inserting this dependence on orbital parameters into our 
previous derivation in Equation ( fTTI i allows us to effectively 
replace the part of 'J'bar which parameterized our ignorance of 
orbital parameters {F(...) in Equation[T9l). giving 

¥b*(0,H,...)cx— ^ -. (26) 

\-rr- cos9 

lid 

In short, our previous derivation applies, but the orbital de- 
pendence is now explicit in 'I'bar- 

Revisiting Figures [7] & [8] recall that we included this or- 
bital dependence in the predicted curves therein. For each 
orbit, the predicted curve is given by the solution for the gas 
mass within the critical Rg^/Rd (Equation \TJl with the de- 
pendence on orbital parameters as in Equation [26]- we insert 
the appropriate orbital inclination 8 and impact parameter b 
for the two disks in the orbit and sum their expected /burst or 
/disk to derive the model prediction. The difference between 
the different orbits does not appear dramatic in Figure|7]- but 
this is because the burst fraction /burst is plotted on a linear 
scale, suppressing the dependence on 8 at small / gas (most of 
the visible dynamic range in the plot is at large / gas - in this 
regime, however, the stellar bar is weak in any case because 
there is not much stellar mass in the disk - so the result is 
that much of the gas survives and becomes part of the disk, 
regardless of orbital parameters). 

However, the difference between different orbits is much 
more clear in Figure [8] displayed on a logarithmic scale. At 
low / gas , there is much more stellar mass in the disk than gas 
mass, so in principle the stellar bar could (if maximal) easily 
torque away all the angular momentum of the gas. Here, how- 
ever, the orbital parameters become important in determining 
just how efficient this process should actually be. For the or- 
bits close to retrograde (cos8 w -1), the scaling we have just 
derived suggests that ^bar should be suppressed by a factor 
~ 2. But for orbits close to coplanar prograde (cos 8 w 1), 
\l/bar is enhanced by a factor ~ 2- 3 - in other words, because 
the orbit is nearly resonant, the effective co-rotation resonance 
(the orbit interior to which the gas can efficiently lose angular 
momentum to the induced stellar bar) is moved out by a sub- 
stantial factor, including a larger fraction of the disk gas (in 
those extremes, only the gas at very large radii survives the 
merger). 

Again, Figures |7][8] demonstrate that the simple scalings 
based on our model provide an accurate description of the be- 
havior in the full numerical experiments. Figure [9] combines 
these into a single plot - we compare the disk fractions in 
our simulations to the full expectation based on our derivation 
thus far as a function of both gas fraction and orbital parame- 
ters. As noted above, the agreement is good, with a reasonably 
small scatter. 

These results are for four representative orbits spanning a 
reasonable range in orbital parameters - those for which we 
have a large number of simulations covering a wide range in 
the space of other parameters. We consider them first because 
this allows us to robustly determine that the predicted orbital 
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FIG. 10. — Effective efficiency of bars (the parameter *J/, efficiency at 
torquing gas into a starburst, removing its angular momentum and destroying 
the disk), as a function of the effective orbital parameter. Each point repre- 
sents the effe ctive constraint on 9? from fitting a correlation of the form in 
Equation )19t to a suite of simulations over a range of gas fractions, masses, 
and mass ratios (since we are interested in a comparison of orbital parameters 
here and not e.g. gas fractions, we choose to normalize so that < \& < 1). 
Black points are the most well-sampled orbits (e, h, k, f), shown in Figure[8] 
purple points are more limited studies of orbits (c, i, m, d), red points are 
a study of major and minor mergers with 4> = 0, 9 = 0, 30, 90, 150, 180°. 
Top: <E" c ff versus mean orbital inclination (cos (9)) (for 1:1 mergers, we av- 
erage the two inclinations). Black solid line is the simple linear scaling 
V^eff °c 1/(1 — 0.6(cos#)) from Equation (26); red dot-dashed lines are the 
numerical solutions for the appropriate (cos6>), as in Figure[8] Middle: ty e ff 
versus the full numerical expectation (properly solving for the effective bar 
strength in each disk and then adding the burst fractions, rather than just tak- 
ing \I'((cos9)). The efficiency of disk destruction and angular momentum 
loss scal es with orbital parameters in the simple manner predicted in Equa- 
tion (26). Over a typical random cosmological ensemble of orbits, we expect 
values similar to those between our typical e and f orbits. Bottom: Vf^ versus 
net specific angular momentum of the system (adding/cancelling the initial 
disk plus orbital angular momenta and dividing by the final baryonic mass). 
The two are actually anrf-correlated, demonstrating that disks do not arise 
after or survive mergers owing to co-addition of angular momentum (and 
cancellation is inefficient at destroying disks) - rather, systems with aligned 
angular momentum vectors are in greater resonance, triggering stronger in- 
ternal asymmetries in the primary that drain more angular momentum from 
the gas. 



scalings do not depend on e.g. stellar mass, halo properties, 
feedback prescriptions, or other varied physics in the simula- 
tions. Having done so, we consider a more limited sampling 
of a much broader range in orbits given by TableQ]in order to 
survey the full dynamic range of orbital parameters. 

Figure [10] shows the results of this. For a given suite of 
simulations with some particular orbital parameters, we first 
construct the correlation /disk(/gas) as in Figure[8] Rather than 



adopt some a priori model for the orbital dependence, we then 
fit the points in that correlation to a function of the form in 
Equation ( [T9T > - i.e. effectively fit for the normalization or "ef- 
ficiency" of angular momentum removal, which we define as 
(^eff(orbit)). We compare this, for our ensemble of orbits, 
to our analytic expectation from the simple scaling in Equa- 
tion (l26l l and to a full numerical solution (technically, for 1 : 1 
mergers, we want to solve this separately for each disk and 
add the two, although just considering the primary is a good 
approximation for less major mergers). The agreement with 
our analytic model is quite good across the entire range of 
orbital parameters, implying that we have captured the most 
important physics of resonant interactions in this simple scal- 
ing. 

We also compare this effective efficiency of disk destruc- 
tion with the net specific angular momentum of the merger 
remnant (assuming pure addition/cancellation of the initial 
baryonic angular momenta of the disks and the orbital angular 
momentum). The result is actually an anti-correlation: sys- 
tems with aligned angular momentum vectors, e.g. coplanar 
prograde mergers being the extreme case, induce the most ef- 
ficient bars and remove angular momentum most efficiently 
from the gas. Systems where the angular momentum vectors 
are misaligned (e.g. polar orbits) or anti-aligned (retrograde) 
actually leave the largest disks in place. This clearly empha- 
sizes that it is not, in fact, any direct addition/cancellation of 
angular momentum that determines or enables disks to form 
in and survive mergers. Rather, the cases with the largest net 
angular momentum are most resonant, inducing the strongest 
resonant asymmetries in the merging pair, which most effec- 
tively drains angular momentum from the gas and leaves a 
compact, bulge-dominated remnant. 

4.3.3. Dependence on Mass Ratio 

The major remaining parameter to study is the merger mass 
ratio. Thus far, we have restricted our attention to equal mass 
1 : 1 mergers, which allowed us to make several convenient 
simplifying assumptions. Nevertheless, most of our previous 
derivation applies. None of the scalings that we have explic- 
itly derived up to now are dependent upon mass ratio. How- 
ever, we have quantified the strength of the induced stellar 
and gas bars with the parameter ^bar, which we expect should 
scale with mass ratio. Moreover, we have made the assump- 
tion that the pre-merger disk stars are entirely violently re- 
laxed by the merger. While this is a good assumption for 1 : 1 
mergers, it is not true for minor mergers (a 1 : 10 mass ratio 
merger, even with no gas, will clearly not transform the entire 
primary stellar disk into bulge). 

First, consider this stellar component: there are a number of 
ways to derive the disturbance of the stellar component in the 
merger. The simple expectation is that the mass in galaxy M\ 
which can be violently relaxed by collision with galaxy M2 is 
proportional to M2/M1 - the net energy deposit, tidal forces, 
and the mass fraction brought in from a potentially disrupted 
satellite all scale in this manner. For simplicity, consider the 
case where the secondary M2 is much smaller and more dense 
than the primary, and falls in on a nearly radial orbit in the 
final encounter (which is a good approximation, given the ef- 
ficiency of angular momentum transfer from the orbit to the 
halo). Since we are assuming M2 <C Mi, treat M2 as a point 
mass, and consider its final orbital decay, where it oscillates 
with rapidly decaying amplitude through the center of the pri- 
mary with initial impact velocity v; 1=3 v c and damping spatial 
amplitude £ max < R c /. At some instant, then, the secondary is 
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at location (R',<f)', z') = £(co&9, 0, sin#) (we rotate such that 
the secondary orbit defines = without loss of generality). 
A star in the primary disk at (R,(j>,z = 0) then feels some po- 
tential from the secondary (= $2) and experiences a vertical 
deflection out of the disk d® 2 /dz = (GM 2 /R y )£ sm6f{l/R), 
where f(u) = [1 +u 2 -2ucos9 cos0] 3 / 2 ~ 1. We are only in- 
terested in the time the secondary spends at I ~ R (when its 
much closer to the disk or further away, the vertical pertur- 
bation is weak), so it effectively acts for a time St ~ R/vi as 
it passes through I ~ R in its ringing about the center. If we 
know the full potential, we can solve for the deflection as a 
function of time and calculate the full acceleration of the disk 
stars at (R,<fi), which yields an effective net velocity deflec- 
tion while the secondary is on one side of the galaxy of Sv = 
(GM 2 / Rvi)(£ max s'm9 / R) f(£ m:ix / R). Deflection occurs when 
Sv ~ v or larger, so if v = v(r) v c (where v depends weakly 
on r) and we substitute for v c = *JGM\ /Rd here and in Vj we 
obtain the criterion GM 2 /R > v c v, ~ v 2 c ~ GM\ /Rj, i.e. (rear- 
ranging) R/Rd <M 2 /Mx sin6> = fj, sind. The sin# dependence 
comes because we considered only vertical deflection of stars 
(i.e. some heating to v 2 )- a coplanar orbit (in this limit) will 
obviously induce no such heating, but will introduce deflec- 
tions in the radial direction (heating v|). We can repeat our 
derivation considering where these deflections are significant, 
and find (as one would expect) R/Rd < [i cos 9. So, the ab- 
solute mass fraction scattered should be more or less angle- 
dependent, although the orbital anisotropy (3 Z = \ -v\ jv\ will 
depend significantly on the orb i tal inc lination 9. 

For the case of a thin M estell (1 19631) disk with no bulge, we 
can solve these equations exactly and obtain the simple so- 
lution that a merger with a secondary M 2 scatters exactly M 2 
worth of stars in the primary, completely independent of the 

inclination 9, (but with an anisotropy (3 Z (9) ~ 1 - Tf^^ej - 
although this ignores a proper treatment of further mixing as 
the perturbed stars interact with each other, and thus does not 
reproduce orbits quite as radial as seen in simulations). The 
full numerical solutions for arbitrary cases yield the general 
result that, when inside a radius that encloses a mass ~ M 2 
in the primary, then the presence of the mass M 2 is a signif- 
icant perturbation, which scatters those stars in the primary 
- i.e. deflections occur rapidly, so the stars violently relax. 
At larger radii, where M enc — * M\ > M 2 , the motion of the 
M 2 secondary at the center is a small perturbation. The disk 
at these radii is perturbed adiabatically by the motion of the 
secondary, which can induce some warps and/or disk heating, 
but will not violently relax the stars. Reversing this derivation 
for the secondary, it is trivial that essentially all the mass in 
the secondary (we ignore stripping of the tightly bound stellar 
mass) will be violently relaxed. So the total mass violently re- 
laxed will be ~ M 2 (in the primary) plus M 2 (the secondary), 
out of a total mass Mi +M 2 - i.e. in terms of the mass ratio 
/i = M 2 /M\, the fraction of the pre-merger stellar disk mass 
which is destroyed and turned into bulge is 

2 11 

/*.disk(destroyed) = . (27) 

1 +/i 

Technically this assumes the systems are initially pure disk, 
but the corrections if they have pre-existing bulges are not 
large (generally smaller than the simulation-to-simulation 
variation; although we discuss them in more detail in § |5), 
so this is a reasonable approximation for general cases. 

Now, consider the gas. It turns out that a similar linear scal- 
ing in Equation d27T i is found for how the gas mass in the 



starburst (i.e. the fraction that loses its angular momentum) 
scales with mass ratio, as one might expect. In detail, though, 
the derivation must be revisited (and will include additional 
terms depending on orbital parameters): because the gas is 
collisional, even a large vertical deflection of gas at some R 
does not translate to a loss of that gas disk, since the gas can 
dissipate the vertical energy and no loss of rotational angular 
momentum has occurred. Deflections in the R direction will 
be resisted by hydrodynamic forces. So, for a proper deriva- 
tion, we return to our model of the stellar bar torquing the gas 
bar. The essential question is how the amplitude of the in- 
duced stellar bar (our term ^bai) should scale with mass ratio. 

Take the thin disk limit (this is just for convenience, the 
final scaling is unchanged if we allow for a finite stellar disk 
thickness); the disk surface density is linear in the potential 
according to Poisson's equation, 

V 2 $ = 47r£(fl,0)c5(z). (28) 

So, since the non-axisymmetric potential of the secondary, at 
some distance b (the impact parameter), must scale as roughly 
$ ~ GM 2 /b 3 , we expect the amplitude of the induced bar 
(perturbation in E) should also scale as M 2 /b 3 . Fractionally, 
this yields ^ b ar «Mb ar M ~ fi(h/b) 3 . 

We can show this more formally: if <&o is the (azimuthally 
symmetric) potential of the primary and $1 is the pertur- 
bative potential of the secondary, which induces the sur- 
face density perturbation Si oc / bar oc \&bar that defines the 
bar, we have V 2 <E>i = 4ttSi(/?, <j))S(z). We can expand any 
potential $1 as & a (kR) exp[i(mcf>-ujt)-k\z\], which gives 
the trivial solution Sj = Yl a (kR) exp[i(mcf)-uit)-k\z\] where 
S fl =-\k\/(2nG)d 2 <i> a /dr 2 . We expect <S> a - -GM 2 /r, so we 
obtain d 2 & a /dr 2 ~ -2GM 2 /r i (note that we can generalize 
this to extended distributions for the secondary, relevant for 
e.g. more major mergers, and the change is, for reasonable 
profiles, equivalent to replacing r with Vr 2 + a 2 ). The de- 
tails of the mode structure turn out not to be important, since 
the behavior we are interested in is dominated by modes with 
\k\ ~ l/Ra; but we can, for example, treat Q a (kR) as the po- 
tential generated by a point source (appropriate for e.g. the 
small mass-ratio limit) at the impact parameter b and expand 
the wave modes appropriately, then integrate over the modes 
in the disk to determine the bar strength (i.e. the total mass ef- 
fectively contributing to the bar). In any case, up to a numeri- 
cal constant that is weakly sensitive to the mode structure, we 

obtain M bar oc (fc2+ ^ )3/2 M 2 . In terms of * bar = M bar /M u this 
gives 

* bar ocM(l + [W 2 r 3/2 , (29) 

where again b is the impact parameter and the 1 + [b/ R^] 2 term 
effectively allows for the interpolation between the case of 
an orbiting point mass an d a penetrating encounter (see e.g. 
iBinnev & Tremaindll987l) . 

As noted in § 14.3.21 if we are just interested in the end prod- 
uct of a merger - i.e. we do not care what happens on each pas- 
sage separately as the companions lose angular momentum, 
but only in the surviving disk fraction and total burst fraction 
- then we are not interested in some initial impact parameter 
b but only in the impact parameter on the final passages close 
to coalescence, when angular momentum loss has made the 
orbits nearly radial and the forcing is strong. We can see di- 
rectly from Equation d29b that the forcing is dramatically sup- 
pressed by a factor ~ (b/Rd) 3 on earlier, large-impact parame- 
ter passages, so these can be effectively ignored in calculating 
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the remnant properties (we confirm this is true in a sample of 
simulations with much larger R per i). Eventually, for any sys- 
tems which are destined to merge, angular momentum trans- 
fer yields a nearly radial orbit with b ~ R&, and this is where 
most of the forcing occurs, so the remnant solution is effec- 
tively given by ignoring the b dependence above (technically 
summing over each passage with the appropriate b is possi- 
ble, but in practice we obtain the same result to within the 
simulation-to-simulation scatter by assuming b — > in Equa- 
tion [29}. The final dependence on any initial impact param- 
eter b is therefore weak, so long as the systems are bound to 
merge. The dependence on mass ratio /i, however, is fixed. 

If the mass enclosed is linear in "fbar (the case for e.g. the 
Mestej d 19631) disk and an isothermal sphere, and not a bad 
approximation for the regime of interest for an exponential 
disk), we then have a similar result for the gas as the stellar 
distribution: the secondary induces a burst of mass oc M2 in 
the primary M\ . Reversing the derivation, the primary (since 
it is larger, so M\ /Mi > 1 induces a burst (assuming the two 
have similar initial gas fractions) oc M2 in the secondary (i.e. 
bursting all its gas). The net burst mass oc 2M2 relative to the 
remnant mass Mi +M2 is then 

2 fi 

/bustoc- . (30) 

1 +fi 

This is generally applicable for mergers; however we will note 
below that, because they do not coalesce (and therefore do not 
eventually come in with b — > or brake their orbital energy 
interior to the stellar distribution of the primary), this is not 
exactly applicable to e.g. fly-by or first passage scenarios. 

Figure[T2]tests this prediction in an ensemble of simulations 
spanning a range in mass ratio from fi = 0. 1 - 1 . For a given 
set of orbital parameters (fixed), we plot the disk and burst 
fractions (/disk and /b U i- s t) of the remnant, as a function of the 
immediate pre-merger gas fraction / gas , as in Figures[7][8] For 
the 1 : 1 mergers, we plot our expectation based on the simple 
scaling in Equations ([T7]|-(|26]|, including th e dependenc e on 
/gas and orbital parameters (8) following § l4.3.lE~3~2l We 
then show the prediction for mass ratios 1 : 2, 1 : 4, and 1 : 8, 
according to our derivations here. 

This includes two important corrections: instead of assum- 
ing the entire stellar disk is turned into bulge (which was a 
good approximation for the 1 : 1 mergers), we allow the frac- 
tion of the stellar disk that is destroyed (turned into bulge) to 
depend on mass ratio following Equation d27l i - so some (con- 
siderable) fraction of the disk is assumed to survive in higher 
mass-ratio mergers. We also include the scaling with mass 
ratio in ^bar, used as before to calculate how much of the gas 
participates in the starburst. So, in the high mass-ratio cases, 
both the fraction of the gas that loses its angular momentum 
(fraction of / gas ) and fraction of the pre-merger primary stellar 
disk turned into bulge (fraction of (1— / gas )) are suppressed by 
a factor ~ fj,. 

For each of the orbits surveyed (and the range in e.g. abso- 
lute masses, gas fractions, and feedback prescriptions in our 
minor merger simulations), this simple rescaling according to 
the merger mass ratio provides a good approximation to the 
behavior in the full hydrodynamic experiments. Both the total 
surviving disk fraction (which reflects both the ability of the 
pre-merger stellar disks and the pre-merger gas to survive the 
merger) and the burst fractions (which reflect only how much 
of the gas survives/loses angular momentum) are accurately 
predicted, suggesting that our derivations are reasonable for 
both the dissipational and dissipationless components of the 



galaxy. 

Figure [12] summarizes these results. We first compare the 
final disk fraction in the simulations to our prediction includ- 
ing e.g. the dependence on gas content and orbital parameters 
but without any accounting for mass ratio (assuming all merg- 
ers are just as efficient as a 1 : 1 merger). Unsurprisingly, this 
works for the 1 : 1 mergers, but is a terrible approximation 
to mergers of very different mass ratios. We then compare 
allowing for the same scalings but including the predicted 
mass ratio dependence. The agreement between full simu- 
lation and our simple analytic expectations is good - with a 
scatter for the high disk fractions typical of intermediate and 
minor merger remnants as low as ~ 20%. 

One important caveat here is that, for mergers of increas- 
ingly small mass ratio /x, the merger timescales become 
long. At the smallest mass ratios we consider, ^1:10, this 
timescale may become sufficiently long that the secular (i.e. 
self-amplifying) instability/response of the disk may become 
important over the duration of the merger. It is not entirely 
clear what the response of such an (initially driven) sys- 
tem will be; whether or not, for example, the driven non- 
axisymmetric modes will remain locked to their driver (the 
secondary orbit) or de-couple and move at the pattern speed 
dictated by the internal stability properties of the disk. This 
competition between secular processes (more sensitive to e.g. 
the detailed structure, rotation, and pressure support of the 
disk) and merger-driving in this regime probably contributes 
to some of the increased scatter in burst fractions seen in Fig- 
ure QT| at the lowest /i. For this reason, it is reasonable to re- 
strict a definition of "mergers" to this mass ratio and more ma- 
jor interactions: at smaller mass ratios, secular/internal pro- 
cesses (even if initially driven by interactions) may be more 
important than the direct driving from the interactions them- 
selves (or at least operate on comparable timescales). 

4.3.4. First Passage and Fly -By Encounters 

We have derived a general equation for the disk mass that 
should be lost in mergers, and demonstrated that it is robust to 
variations in a wide range of galaxy properties. Most of our 
derivation is completely generalizable as well to encounters 
where the systems will not merge (or at least are not immedi- 
ately merging). Two cases of interest (which are, in the short 
term, essentially equivalent) arise: first passages and "fly-by" 
encounters (in which there is a close encounter but the veloc- 
ities are sufficiently large to delay or prevent a merger). 

In such a passage, there is of course no violent relaxation 
and mixing of stars, so we assume the stellar disk is left intact 
(excepting the bar response). The same physics will govern 
bar formation and loss of gas angular momentum. The pri- 
mary difference is the suppression by the appropriate impact 
parameter b in Equation d29l . We argued before that the term 
oc [1 + (b/Rd) 2 ]^' 2 should ultimately be neglected for mergers 
because in the final passage(s) that dominate, the condition of 
merging more or less guarantees b — > 0. However, clearly this 
is not the case on a non-merging passage. 

This introduces a non-trivial uncertainty - we quote [1 + 
(b/Rd) 2 ] 3 / 2 where b is the impact parameter and R^ is some 
characteristic scale length of the system. But in detail, the 
appropriate "impact parameter" is really the actual distance 
of closest approach, which is usually somewhat smaller than 
the distance of approach estimated from infinity (the formal 
impact parameter definition), as some angular momentum is 
already lost. Moreover, in detail, is the appropriate R^ the ex- 
ponential scale length? The half-mass radius? Any such radii 
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FIG . 11. — Top: Surviving disk mass fraction as a function of pre-merger gas fraction / gils , for a series of mergers of varying mass ratios (symbols). Each 
panel shows a series of mergers with different orbital parameters (orbits m000-ml80). Color encodes mass ratio of the merger: 1:1 (red), 1:2 (green/cyan), 1:4 
(blue), 1:8 (purple). Dotted lines (of corresponding color) show our prediction (Equation[7) for the given orbital parameters and mass ratio. Note that for minor 
mergers, > f glis is allowed, because some of the original stellar disks are predicted to survive the merger as well as some of the gas which does not lose 
angular momentum. Bottom: Starburst mass fraction in mergers of a given mass ratio p, relative to our model prediction for 1:1 mass ratio mergers with the 
same orbit and pre-merger gas content (symbols, as top panels). Lines show our simple linear model (solid blac k; this does well for typical orbits but the bursts in 
nearly prograde orbits - mOOO - are somewhat more efficient thanpredicted owing to the effects described in § 14.3.6) and full numerical calculation (dot-dashed 
red; two lines correspond to different mass profiles), as in Figure[8| Minor mergers induce less efficient bursts, and do not completely destroy the primary disk: 
the scaling of these efficiencies with mass ratio agrees well with our dynamical model predictions. 
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FIG. 12. — Predicted (as a function of orbital parameters, pre-merger gas 
content, and mass ratio) and actual post-merger surviving disk fraction for 
the simulations in Figure 011 (symbol type and color denote orbit and mass 
ratio in the same style). Top: Comparison assuming there is no dependence 
on mass ratio (i.e. treating all cases as p = 1). Clearly, this is inappropriate for 
minor mergers, but it is also inappropriate for even intermediate major merg- 
ers (mass ratios p = 0.3-0.5). Bottom: Comparison including the predicted 
dependence on mass ratio of both destruction of the stellar disk and angular 
momentum loss in the gas. Our predictions as a function of mass ratio, orbital 
parameters, and gas fraction are accurate in the simulations to ~ 0.1 in /disk- 



are of course closely related, and all of these uncertainties in 
the exact definition change the term b/Rd only at the factor 
~ 2 level, but since the dependence ~ (b/R$f is fairly strong, 
this is important on a quantitative level for these fly-by situa- 
tions. 

In practice, we find that using the impact parameter b de- 
fined as the halos approach (i.e. neglecting detailed resonant 
loss of angular momentum) and taking to be the half-mass 
radius of the system works well in a mean sense. The re- 
sults of this exercise are shown in Figure Q~3] We plot the 
fraction of the gas available at the time of a first passage or 
fly-by encounter which is consumed in the induced burst (we 
define the strength of the induced burst by integrating the star 
formation excess over the interpolation between the pre- and 
post-flyby star formation rates; see e.g. lCox et al. 20081, for 
details), as a function of the gas content, for different orbits 
as in Figure [8] We predict that the efficiency of channeling 
gas into the burst should scale as ~ (1 -/ gas ), as before, and 
that the scaling with orbital parameters should be similar. We 
also show, for cases with otherwise identical gas fraction at 
the time of first passage and the same orbits as Figure [TU 
how this scales with merger mass ratio (again, expected to 
be the same as that we derived above). Altogether, adopting 
our previous estimates, but re-normalizing appropriately for 
the impact parameter of the passage (here i/fij k 1, as we 
defined it above) yields a good approximation to the typical 
behavior in our simulations. We also test the behavior as a 
function of impact parameter, and find that our simple scal- 
ing is a reasonable approximation, yielding rapidly diminish- 
ing bursts as the impact parameter is increased to b ^> R^ (at 
some point here, our estimates from the simulations become 
ambiguous, as a ~ 1 % enhancement in star formation is be- 
low the level of random fluctuations in isolated disks). We 
have also checked whether or not the pre-flyby stellar disks 
are destroyed - as expected, they are left more or less intact 
by fly-by encounters. The disks may be heated, and in fact 
some "pseudobulge" can form from the buckling of the bar 
induced in the stars, but we are not attempting to predict or 
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FIG. 13. — Comparison of the burst fractions in single fly-by (or equivalently, first-passage) scenarios and the general application of our model scalings. 
In these cases, there is no significant violent relaxation (no stars merge), so the stellar disk is left completely intact. Some pseudobulge may result from the 
induced bar and disk heating, but we are not modeling that here. Some burst still results from the same induced non-axisymmetry in the primary, which should 
be described by our same scaling (Equation|7J. The important differences from a case that will merge are: (1) the suppression of the induced burst by a factor 
~ [1 +(b/Ri) 2 ] 3 / 2 (where b is the impact parameter and we find decent agreement with our simulations when is the half-mass disk radius), whereas in cases 
that will merge b — > is appropriate, (2) the lack of violent relaxation of the stellar disk, and (3) an expected increased scatter, as the details of the approach are 
more important (and there is no merger/in-spiral, which tends to average out the exact details of the approach). Top: Burst fraction (relative to gas supply at the 
time of the passage) versus gas fraction. Our simple linear model prediction (black solid) and numerical predictions (red dot-dashed) are shown, with the results 
from the first passages and fly-by encounters of the simulations in Figure [8] appropriate for each set of orbital parameters shown. These cases had b/Rd rj 1 
(and that was used in the predictions - the curves assuming b = 0, as we used for the post-merger systems, would be a factor ~ 3 higher, in conflict with the 
simulations). Middle: Same, but as a function of mass ratio for systems in Figure fTTI with otherwise equal orbital parameters and gas fractions at the time of 
passage. Bottom: Same, but as a function of impact parameter for 1:1 mergers with /g as ~ 0.2 (black) and 0.4 (blue). Note that "burst" fractions < 1% of / gas 
are essentially equivalent to zero (equivalent to random fluctuations in isolated disks). Our predictions describe first passages and fly-by encounters reasonably 
well, although there is larger scatter about them owing to differences in the details of how the passage proceeds. 



study pseudobulge formation here (rather considering it, as is 
often the case in observations, to be fundamentally still part 
of the stellar disk rather than part of a violently relaxed "clas- 
sical" bulge). In an average sense, then, our derived scaling is 
generally applicable. 

However, the details of exactly how the approach proceeds 
will introduce considerable scatter in the amount of burst trig- 
gered on first passages and in fly-by encounters. This is plain 
in the large (factor ~ a few) scatter in Figure [13] Further, 
details such as the structure of the bulge are increasingly im- 
portant in the limit of weak interactions, where distortions in 
the potential of the primary that would trigger gas inflows can 
be suppressed by the presence of a larger bulge (and note the 
caveat from § 14.3.31 that the secular/internal response of the 
disk will become relatively more important in weaker inter- 
actions with smaller mass ratios and larger impact parame- 
ters). We therefore expect in general that our predictions can 
be quite broadly applied, but are less robust for any specific 
case if it is a single fly-by as opposed to an integration over a 
full merger. Fortunately, in the case of systems that will ac- 
tually merge, these details tend to average out or be unimpor- 
tant, yielding the relatively small scatter we have seen in our 
previous predictions. In those cases, we do not need to be too 
concerned with the exact details of the impact approach, nor 
the structural details of the galaxy (in particular because our 
predictions are for integral quantities at the end of a merger, 



various effects will tend to cancel out - for example retain- 
ing more gas on first passage will yield a larger supply for the 
second burst, etc.). 

4.3.5. Independence from "Feedback" Physics 

Our derivation of the torques causing gas to lose angular 
momentum in mergers is purely dynamical. All else being 
equal (i.e. for systems with the same gas content and dynam- 
ical structure at the time of the final merger), we therefore 
expect that the detailed physics of e.g. "feedback" from su- 
pernovae, stellar winds, and AGN activity should make little 
difference. 

Figure [14] demonstrates that this is indeed the case. We 
compare the starburst and surviving disk gas fractions of 
merger remnants, relative to those predicted by our simple 
dynamical model as a function of the merger mass ratio, or- 
bital parameters, and gas content at the time of the merger, 
for suites of simulations with two different prescriptions for 
supernovae feedback and the effective equation of state of 
the ISM. In terms of our q eos parameter (see § |2), we com- 
pare q eos = 0.25 simulations (a nearly isothermal equation of 
state with effective temperature ~ 10 4 K) t o q eos = 1 simula- 
tions (the "full" stiff Springel & Hernquist (2003) equation of 
state, with effective temperature > 10 s K at the densities of 
interest here). There is no significant systematic offset be- 
tween either the median result or the scatter about our simple 
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FIG. 14. — The effects of feedback on disk survival in mergers. We show 
the distribution in /burst //burst, pred. i e. the burst mass fraction, relative to that 
predicted (or equivalently, the mean in our simulations) for the given pre- 
merger gas fraction and orbital parameters, for simulations with two different 
ISM gas feedback prescriptions (different effective equations of state g C os). 
We also show the corresponding (but measured differently) disk mass frac- 
tions /disk/ /disk, pmd- There is perhaps a small offset in the sense expected (a 
suffer, higher-feedback equation of state for the ISM suppresses bursts by an 
average factor ~ 1.1— 1.2), but this is much smaller than the simulation-to- 
simulation scatter. For a given gas content at the time of the merger, then, 
feedback makes almost no difference (true for AGN feedback and starburst 
winds as well). 
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FIG. 15. — As Figure Ua\ but comparing the distribution of disk fractions 
in merger remnants relative to our simple predictions as a function of stel- 
lar wind and quasar feedback prescriptions. We plot the distribution of disk 
fraction /j; s k in simulations relative to our predicted /diskC/gasjA 4 ,*?) (i.e. our 
calculation as a function of immediate pre-merger gas fraction, orbital pa- 
rameters, and merger mass ratio). Left: Varying starburst-driven wind pre- 
scriptions. We compare our usual weak stellar wind scenario (winds self- 
consistently can generate from the hot gas, but additional mass loading is 
only ~ 1% of the star formation rate) to a fast winds scenario (with addi- 
tional mass loading ~ 0.5M* and wind launch velocity ~ 800kms ) and 
a slow winds scenario (additional mass loading ~ 2M„ and launch velocity 
~ 200 kms~' ). Right: Simulations with and without feedback from accreting 
black holes. For otherwise fixed merger parameters (orbit, mass ratio) and 
disk properties (cold gas content, mass profiles), these feedback prescriptions 
make no different to the starburst or surviving disk gas fractions. They will, 
however, change the the gas content, consumption, and distribution leading 
into the merger. 

analytic expectation. At most, there may be a ~ 20% system- 
atic offset, in the sense that more highly pressurized systems 
(<7eos = 1) have slightly more gas survive - a small offset like 
this is expected because the bars in these cases are slightly 
more "puffy" (essentially the same as a slightly thicker disk 
- for which we derive an analytic expectation in Equation l20l 
that yields an expected ~ 10-20% difference at most based 
on the full possible range of q e0 s). In any case, such an off- 
set is small relative to other systematic uncertainties in disk 
structure and the scatter about the median predictions. 

In Figure Q3] we perform a similar exercise for cases 
with and without central supermassive black holes (we have 
also examined initial BHs with varying initial masses from 
< 10 5 M Q to ~ 10 7 M Q ), and cases with or without a sim- 
ple implementation of starburst-driven winds where winds 
are launched (in addition to the stellar feedback implicit in 
our multi-phase ISM model) with a mass-loading efficiency 
M w i D d = rj w M* relative to the star formation rate and en- 
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ergy loading ef ficiency e w relative to the total energy (for a 
ISalpeterl (1 19551) IMF) available for supern ovae (sampling th e 
range r] w ~ 0.01 - 10 and e w ~ 0.1-1; see lCox et alj |2008b). 
Shown in Figure Q3] are our fiducial weak winds (r) w ~ 0.01, 
e w ~ 0.0025), cases with moderate mass loading into very 
fast winds (rj w = 0.5, e w = 0.25, yielding a wind launch speed 
~ 800 km s" 1 ), and cases with high mass loading but cor- 
respondingly slower wind velocities (r/ H . = 2.0, e w = 0.0625, 
yielding a wind launch speed ~ 200 kms" 1 ). In all these cases 
we find a similar result: at otherwise fixed properties at the 
time of merger, feedback makes no difference to our conclu- 
sions. 

The reasons for th is are described in § 14.3.21 and in more 
detail below (§ 14.3.71 ). Recall, the distortion in the primary is 
driven by the secondary and as such depends only on the grav- 
itational physics of the merger. Given this distortion, the grav- 
itational torques within some characteristic radius are suffi- 
ciently strong to remove the angular momentum from the gas 
in much less than an orbital time. Feedback, then, insofar 
as it changes the effective pressurization or equation of state 
of the gas or drives a wind, is largely irrelevant: because the 
angular momentum is removed in a timescale much shorter 
than the orbital time, the gas (regardless of the strength of 
feedback) cannot dynamically respond with these hydrody- 
namic forces, but must essentially free-fall into the center of 
the galaxy where the starburst is triggered. The radius inte- 
rior to which the torques are strong is not a function of e.g. 
the stability of the galaxy to perturbation, because it is not an 
instability in the first place, but a driven distortion in the sys- 
tem. Moreover, the entire system is strongly in the non-linear 
regime for the time of interest (when we consider interactions 
of the magnitude simulated: mass ratios ^1:8 and more ma- 
jor mergers) - no amount of making the system more robust 
against linear instability would be sufficient to avoid a strong 
gravitational distortion in the violent coalescence surround- 
ing the actual merger (this must be so, because the distortion 
occurs where the disturbances in the potential are greater than 
order unity - for hydrodynamic forces to resist distortion, they 
would have to be stronger than large-scale gravitational forces 
in the equilibrium system, negating the concept of a rotation- 
ally supported thin disk). So what matters is instead where 
that coalescence occurs and how long it introduces such a 
strong distortion, relative to e.g. the local dynamical or orbital 
time of some disk element, giving rise to the simple dynami- 
cal criteria for angular momentum loss developed here. 

That is not to say that for fixed initial conditions (signifi- 
cantly pre-merger or e.g. at first passage), feedback will not 
change the result. There are two primary means by which 
feedback can indirectly have a strong influence on disk sur- 
vival: 

(1) Retaining Gas (Lowering the Stellar Mass Frac- 
tion): As h as bee n demonstrated in a number o f works 
(IWeiletalJ [l99l ISommer-Larsen et ail fl999l 120031: 
iThacker & Couchmanl 120001 I200U iGovernato et alJ 
| 2007t iRobertson et all l2006ab ISpringel et al I l2005b | : 
Springel & Hernquistl 120051: lOkamoto et alJ 120051 
IScannapieco et alJ [2008), these forms of feedback can 
have dramatic implications, in even a short time period, for 
the rates at which cooling of new cold gas from the halo 
and consumption of existing gas by star formation proceed. 
In cases with no feedback, star formation may exhaust 
gas efficiently, leading to predicted systems that are much 
more gas-poor at the interesting time of the final merger - 
according to our model, then, these will not be able to form 
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disks as efficiently as more gas-rich systems. In cases with 
strong feedback from e.g. star formation to lower the effective 
star formation efficiency and recycle gas, the predicted gas 
fractions at the time of merger (from some gas-rich initial 
conditions) could be much higher. Inclusion of stellar and 
supernovae feedback responsible for injecting energy and 
turbulent pressure into the ISM may also be necessary to pre- 
vent the onset of clumping and disk fragmentation in isolated 
gas-rich cases, enabling the stabl e existence and evolution of 
quiescent gas-rich disk s (see e.g. lSpringel & Hernquistll2003l : 
iRobertson et al.l 120041) . In short, feedback may be critical 
to give rise to high gas fractions in the first place, which 
we have shown have dramatic implications for the survival 
of disks - but for a given gas fraction (however that comes 
about in the first place), the results of the merger will be (in 
the short term) independent of feedback. 

(2) Changing the Spatial Distribution of Gas ("Kicking 
Gas Out" of R max ): Recall, our derivations demonstrate that 
it is not necessarily a fixed fraction of gas that loses its an- 
gular momentum: rather (see Equation [TTI and § 14.3. Ill4.3.2b 
it is the mass inside some radius /? max /^d relative to that of 
the stellar disk (characteristic radius /?d) which will lose its 
angular momentum. If some form of feedback can change 
the spatial gas distribution, then, it could have dramatic im- 
plications for disk survival. We have used the radius R m!lx to 
estimate the mass fraction that will burst by assuming the gas 
density profile is broadly similar to that of the stars (which 
is true in our simulations, given their feedback prescriptions). 
But one could easily imagine the extreme limit, where some 
strong feedback keeps all the gas at large radii r ^> R max ~ 
(i.e. a case in which there is a large hole in the gas distri- 
bution, or in which the gas is at least much more extended 
than the stellar distribution) - the stellar disk torques only act 
effectively within /? max , so only a tiny fraction of the gas in 
such a case would lose its angular momentum. Especially 
at high redshift, this may be important in avoiding overcool- 
ing and the fo rmation of too much bulge mass in many sys- 
tems (see e.g. IRobertson et al.l l2004t iGoyernato et al.l 120071: 
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2007)^ A gain, we stress that for a given gas density profile at 
the time of merger, our calculations are independent of feed- 
back; but if feedback alters the gas profile - keeping the gas 
at radii ^> /? max , then it will largely survive the merger. 

4.3.6. Exceptions and Pathological Cases 

We have derived a general model for how disks are de- 
stroyed in mergers and shown that it applies to a wide range 
of gas fractions, orbital parameters, galaxy mass ratios, and 
prescriptions for feedback and gas physics. However, there 
are some pathological cases of more than academic interest, 
as these can explain some small differences with previous re- 
sults as well as illustrate the important physics in our model. 

For example, consider the starburst mass fraction and sur- 
viving disks in our "h" orbits: i.e. a prograde-prograde, copla- 
nar merger of two disks. In this case, the angular momentum 
vectors of both disks and the orbital angular momentum are all 
perfectly aligned. Naively, one might then expect that these 
unique cases would create the largest disks. In fact, the op- 
posite is true. This is largely for the reasons we outline in 
§ 14.3.21 - the alignment of angular momentum vectors means 
that the system is in near-perfect resonance, so it excites the 
largest tidal and bar asymmetries that rapidly drain the gas of 
all angular momentum. As we have shown, the much larger 
space of less-aligned orbits is in fact more favorable to disk 



survival. 

However, while the perfect resonance means that the bar 
efficiency "f^r is large, the amount of mass in the stellar bar 
still scales as 1— / gas , so the burst fraction should vary as 
/gas(l ~/gas) in our simple model (i.e. we would still expect 
that a 100% gas disk would have no stellar bar, hence no burst, 
as we have seen for more representative orbits in Figure|7]i. In 
fact, though, we typically find in these cases that the burst 
fraction seems to scale as /burst = /gas all the way to high val- 
ues of / gas - in short, almost all the gas always bursts - there 
is no suppression by a 1 — / gas factor as would be expected if 
the stellar bar were doing the torquing. 

The reason for this is simple - again, the orbits here are 
perfectly coplanar and in resonance; so this is the one case 
where the secondary galaxy as a whole can directly act as an 
efficient torque on the gas. In short, because the systems are 
perfectly coplanar and in a resonant orbit, the entire secondary 
galaxy (all baryons and dark matter within the stellar R e ) acts 
directly to introduce a non-axisymmetric potential perturba- 
tion (the secondary itself plays the role of the bar). So be- 
cause of this, to an even greater extreme than our scalings for 
more general orbits would predict, this narrow range of orbits 
is pathological and biased against disk formation. However, 
understanding why this is the case, we can check and explic- 
itly show that it is not so for more general orbits, even nearly 
prograde-prograde orbits (such as case e) - in all those cases, 
even those just slightly out of coplanar resonance, the stellar 
bar is indeed the primary source of torque, and our assump- 
tions are justified. This example therefore nicely illustrates 
what the consequences would be if our fundamental assump- 
tions were not true, as well as showing why they are in fact 
true for non-pathological cases. 

Another pathological case of interest is one in which the 
disks are 100% gas at the time of merger. Here, as we have 
said, our simple model predicts no starburst or angular mo- 
mentum loss. In practice, there will still be some loss of an- 
gular momentum owing to direct cancellation in e.g. shocks 
between the disks; but as discussed in § 14.21 there will also be 
the possibility of some gain owing to the angular momentum 
of the merger. In fact, over the range in mass ratios p ~ 0.1 — 1, 
for a range of typical impact parameters b ~ 0.5-5, the ex- 
pected final specific angular momentum from after cancella- 
tion is approximately equal to the initial specific angular mo- 
mentum of the primary (with ~ 20% scatter). Cancellation 
is therefore inefficient. A random distribution of orbits might 
negate ~ 20% of the angular momentum in ~ half the sys- 
tems merging, but will leave ^ 80- 100% of the disk intact. 
Even these cancellations, we find in detail, do not generally 
yield a starburst in the same manner as a merger-induced bar, 
but lead to moderate disk contraction (and an equal number 
of mergers will scatter towards the opposite sense leading to 
disk expansion, keeping a mean specific angular momentum 
that is constant). They do not cause a starburst because, if 
two random parcels or streams of gas shock and lose angular 
momentum, the alignment and relative momenta would have 
to be near-perfect for them to lose, say 95% of the angular 
momentum and fall all the way to the central ~ lOOpc where 
a nuclear starburst would occur. Rather, they will lose some 
fraction of order unity of their angular momentum, fall in to a 
radius smaller by a factor ^2-3 (but not to very small radii), 
and continue to orbit. Without the bar that can continuously 
drain angular momentum, t he true burst is indeed inefficient. 

Although we show in § 14.3.51 that the physics of interest 
are generally independent of feedback prescriptions, there are 
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some patho logical feedback r egimes. These are discussed 
in detail in ICox et all ((2008b); here, we outline the patho- 
logical behavior. If e.g. starburst-driven winds are imple- 
mented with extreme efficiencies M w ; n d 3> M* and with mod- 
erate to large velocities > 200km s" 1 , then there is no defin- 
able "starburst" in the simulations any more, even when the 
gas loses angular momentum - indeed, it becomes almost 
impossible to trigger starbursts by any mechanism. This is 
because the feedback is so extreme that any parcel of gas 
that begins forming stars above some threshold rate is im- 
mediately blown apart and drives away all the surrounding 
gas. However, observations suggest that these cases are al- 
most certainly not relevant - observationally inferred mass- 
loading factors of winds are well below the predicted thresh- 
old where we see this behavior (see e . g. [Veilleux et al.l|2005l : 
iMartinl [1991 120061: lErbelaT] 120061 ISatoe t al. 2008J), and 
moreover the ubiquity of starbursts and recent st arburst rem- 
nants in observed gas-rich major mergers (e.g. ISoiferet al J 
Il984albl: IScoville et all 1 19861 ISargent et alj|1987l 1 1989b im- 
plies that feedback, while still potentially efficient, is not able 
to "self-terminate" a starburst before it even begins (this is in 
fact directly confirmed in observations of out flows in ongo ing 
massive, merger-induced starbursts; see e.g. Martin 2005). A 
similar pathology can appear if we include extreme coupling 
of black hole feedback to the galaxy gas (e.g. allowing 100% 
of the BH accretion energy to couple efficiently), but this is 
also ruled out observationally, both for the arguments above 
(starbursts exist, and the winds seen are not so enormou s; see 
the di scussion in lCox et al.l200 8b; Hopkins et al. 2OO5<i303 
2006), and because such a prescription yields black hole 
masses orders-of-magnitude discrepant from t he observed 
dFerrarese & MemtflfoOOt iGebhardt et alj|2000l) M BH - o re- 
lation (see e.g. Hopk ins et al . 2007© 

4.3.7. Longer-Lived Perturbations: Relation to Secular Evolution 

Thus far, we have focused on activity during the merger, 
roughly defined as the short timescale <~ 10 8 yr following first 
passage and coalescence. In this regime, we have shown that 
(for typical conditions), the dominant source of angular mo- 
mentum loss is the torque on gas from stars in the same disk. 
However, it i s well known from studies of isolated barred 
galaxies (e.g.lWeinberdll985HHernquist & Weinberd ll 992b; 
Friedli & Bendl993tlFriedli et alJl994lAthanassoulal 2002a; 
Athanassoula & Misiriotisl 120021: IWeinberg & Katzl 120071: 
Kaufmann etail 120071: iFoyle et al.1 l2008h that a long-lived 
bar (regardless of whether the bar is purely stellar or 
purely gaseous) will exchange angular momentum with 
itself (or e.g. gas/stars further out in the disk) and 
the dark matter halo, allowing for further angular mo- 
mentu m loss and building a centr a l bulge or "pseudo"- 
bulge (Patsis & Athanassoula 2000; Athanassoula 2002b; 
iMayer & Wadsley 2004; Berentzen et all2007h . Here, we dis- 
cuss the relation of this process to what we have described 
in our merger simulations: in general, we find that it (while 
potentially very important for the long-term evolution of the 
disk and bulge masses and structure) is a second-order effect 
within the merger itself, and on longer timescales is more ap- 
propriately considered an independent, secular evolution pro- 
cess (despite being initially triggered by a merger), whose 
study is better described in simulations of idealized and long- 
lived bars. 

As discussed in § 14.3.21 there is a limit to how far the 
analogy to barred galaxies can be drawn. Recall, we use 
the term "bar" more generally to represent a quadrupole mo- 



ment or non-axisymmetric distortion in the stellar disk: it does 
not necessarily (and, especially after second passage, usually 
does not) morphologically resemble isolated barred spirals 
and may not even have an m = 2 mode structure. Critically, 
the distortion is driven externally by the gravitational pertur- 
bation of the secondary orbit - it is not th e resul t of an insta- 
bility within the primary. As we note in § 14.3.21 this already 
gives rise to a couple of important distinctions: because the 
distortion is driven by the orbital motion of the secondary, it 
has a characteristic frequency (and corresponding radius) in- 
ternal to which angular momentum loss is very efficient (de- 
termined entirely by the gravitational properties and relative 
motions of the systems, not the subtleties of their internal or- 
bital structure), regardless of properties of the primary (e.g. 
gas phase structure, feedback, etc.) that might otherwise make 
the system more or less stable to the development of internal 
instabilities. 

It is straightforward to estimate the relative importance 
(over a short timescale after the merger) of angular momen- 
tum loss from the gas to the shared stellar bar/distortion in- 
duced by the merger, versus that to itself and the dark mat- 
ter halo (the standard secular scenario, which other than the 
initial driving in the merger, will not be driven by the rela- 
tive gravitational motions but by the more standard bar sta- 
bility and spin-down criteria). Approximating the gas as a 
rigid, thin bar of mass Mb ar , gas ~ /gasMrar and radius R\, m ~ 
Rd, we can estimate the specific torque from the remain- 
ing ga s disk and halo i n the dynamical friction limit, fol- 
lowing Weinberg (1985) (for more detai led solutions, which 
ultima t ely give similar results, see iHernauist & Weinberg! 
1 1992bt I Athanassoulall2003t IWeinberg & Katzll2007l) : d//df = 

-47raG 2 Mbar,gasP(^/)Vbar' wnere P(R(i) is the background 
density and a ~ 1 is a numerical constant (depending 
on the exact shape of the bar, potential, and phase-space 
di stribution of t he background). If the "background" is 
a iMest el (1963) disk or isothermal sphere, this becomes 
-2aGMbar. g as/^rf- Compare this to our Equation ( fT2b for 
the instantaneous torque on the gas bar from the stellar bar: 

-GMbar. */(Rd V sin 2 4>b+H 2 ). Removing the common fac- 
tors, the gas/halo torque goes as ~ / gas , whereas that from the 

stellar bar goes as - (1 - / gas )/ V 'sin 2 <j> b +H 2 ~ (1 -f m )/<Pb 
(because <fii, ~ H <C 1). In short, the torque from the gas disk 
and halo goes as / gas because it is a second-order resonance 
effect (amplified and trading off with the gas bar), whereas the 
stellar bar strength goes as the stellar mass fraction (1 -/ g as), 
but boosted by a factor ~ 1 /4>i, representing the small angle 
of offset between the two bars - i.e. the stellar bar is in much 
closer spatial proximity (in particular in spatial alignment in 
the disk plane) to the gas bar. 

This simple comparison gives a reasonable quantitative pre- 
diction of the relative torques exerted by the halo and stellar 
disk in our simulations. Essentially, we have just re-derived 
the well-known fact that the timescale for a bar to damp its 
own angular momentum via resonant interactions with itself 
and/or the halo is some number ( ~ a few) bar rotational pe- 
riods (|Athanassouia~ & Misiriotis 2002; Athanas soula! 120031 
IWeinberg & Katzll2007l: iKaufmann et al.ll2007f) (each bar ro- 
tational period being ~ 1 — 2 times the disk rotational pe- 
riod), whereas in the typical mergers the gas is drained of 
angular momentum by the much stronger local torques on a 
timescale much shorter than an orbital time, allowing it to 
more or less free-fall into the galactic center. Comparing these 
timescales gives a similar ratio of torque strengths. Obviously, 
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as /gas — ► 1 . the torque from the halo must eventually domi- 
nate, but this will not happen until /* = (1 —/gas) ^ 4>b ~ 0-1 
(given typical bar lags of ~ a few degrees or the ratio of the 
timescales above). 

In practice, such a situation is somewhat contrived (it is 
very difficult to maintain a disk with a true > 90% gas frac- 
tion), and unlikely to be of broad cosmological relevance (we 
have no simulations in this regime with which to compare, in 
fact, because even initially 100% gas disks with low star for- 
mation efficiencies will be < 80% gas by the time of the actual 
merger). However, our Equation (|7]l can be trivially modified 
to include these effects: the (1 — / gas ) term should be replaced 
with a more appropriate (1— / gas + e/i), where e/, ~ <j>b ~ 0.1 
represents the contribution of angular momentum loss to the 
halo and outer gas disk during the merger. This exchange 
of angular momentum, therefore, sets some minimum bulge 
mass (with mass fraction ~ 10%) that would form even in a 
pure-gaseous disk merger. 

By comparing the relative instantaneous amplitude of the 
torques from the halo/gas and stellar disk, we are compar- 
ing how important they each are in the loss of angular mo- 
mentum from the gas over the same (relatively short) merger 
timescale. More important is the fact that, in a gas-rich case 
where the distortion to the stellar distribution may be an inef- 
ficient torque, the gas bar could be long-lived and continue to 
lose angular momentum over longer timescales. It is not nec- 
essarily clear that this would happen, however - a number of 
studies suggest that gas and stellar bars become self-damping 
once a central mass concentration (i.e. a nuclear starburst trig- 
gered by gas inflows, in this ca se) is in place with a mass 
fraction larger than a few per cent (Bournaud & Combes 2002; 
iBerentzen et alj|2003l |20M l2007t lAthanassoula et al.|[2005l) 
(but see also Kaufm ann et al.l l2007). In such a case, we again 
arrive at the conclusion that these processes set a minimum 
bulge mass from a large bar-inducing perturbation, but do not 
dominate the creation of much larger bulges in mergers. 

Regardless of this effect, it is not clear that a bar can sur- 
vive a substantial merger: recall, the distortion following sec- 
ond passage and coalescence resembles a bar only in that it 
introduces a rotating quadrupole distortion in the disk poten- 
tial (allowing us to describe it as a "bar" for analytic conve- 
nience), not necessarily in its structure or longevity (it does 
not necessarily share the orbital "pileup" that allows a bar to 
survive), and moreover the actual coalescence of the galac- 
tic nuclei will disturb any bar structure that may be present. 
Quantitatively, we find our remnants rarely have significant 
long-lived m = 2 modes in the stars or gas - the mode ampli- 
tude tends to damp after merger on a timescale < 10 yr (i.e. 
the free-fall or dynamical time, much slower than the typi- 
cal significant number of orbital times for standard bar self- 
braking). This is similar to the c onclu sions in the bar stud- 
ies of e.g. [Bournaud & Combes (2002) and Ber entzen et al.l 
d2007l) . who find that the combination of the formation of a 
bulge/small central mass concentration from gas inflows and 
the disturbance/heating to the bar itself in interactions pre- 
vents even gas-rich systems from maintaining or very rapidly 
re-forming a bar after a significant merger (as opposed to a 
fly-by passag e, which may more effi ciently induce long-lived 
bars; see e.g. IBerentzen et aTll2004l) . That is not to say a bar 
may not form in the re-formed remnant disk, but such a bar 
would arise in a standard secular fashion, and should be con- 
sidered in the context of the long-term secular evolution of the 
merger remnant. 

If, however, the potential distortion survives the merger to 



form a stable bar, it can certainly be important to the long- 
term evolution of the system and buildup of the bulge. How- 
ever, this case is outside the scope of this paper, and should be 
more appropriately considered as subsequent evolution of the 
remnant (albeit with an initially merger-induced bar). This is 
because the timescale for the bar to lose angular momentum 
and contract is some number of rotational periods - so the 
gas losing angular momentum will slowly spiral inwards in 
some number of orbital periods (turning into stars and possi- 
bly being ejected by feedback as it does so), rather than free- 
falling into a central burst in a time much less than an orbital 
period. The end r esult of such angular momentum loss can 
resem ble a bulge dMayer & Wadsleyl l2004t iDebattista et al.1 
2004), although the expectation of rotational support and 
"disky-ness" in the material lead to it more likely being a 
"pseudo-bulge" typical of secula r processes (ICombes et al.1 
[19901 [ Kuiiken & Merrifield|[l99l lO'NeiU & Dubinsk|200l 
iKormendv & Kennicuttl 12004; Atha nassoulall2005l) . Depend- 
ing on e.g. details of the equation of state, feedback, and 
rotational support of t he gas disk, it may a lso amount to 
steady disk contraction (Debattista et al. 200 6|) o r emergence 
of a t wo-component disk dKaufmann et al.l l2007; Fovl e et all 
2008). A number of effects will be important in this regime, 
including the effects of feedback in pressurizing the disk and 
smoothing out substructure, and the role of accretion and 
mergers in rebuilding the disk as such evolution continues 
(since it is occurring on timescales ~ several Gyr, compa- 
rable to the characteristic timescales for new accretion and 
mergers). 

These effects make it difficult to predict the net effect of 
such evolution. For example, if in a pure gas merger of 
mass ratio /i the specific angular momentum (on average) 
is increased (by addition of specific angular momenta plus 
orbital angular momentum) by an amount ~ e m n jdisk, but 
then the induced bar (of amplitude ~ /i) loses its angular 
momentum (~ /i ./disk) on a timescale ~ N t mt , the sequence 
of mergers and induced bars compete (given the cosmologi- 
cally ex pected timescale ~ jitu between mergers of mass ra- 
tio \x in iFakhouri & Mai d2008l) . and that for a disk of mass 
fraction relative to the halo, f rot ~ nidtu, one obtains 
dj/dt ~ j ^di s k / 'tn [e,„ - fi/Nma] - i.e. more major events will 
tend to lead to angular momentum loss in gas, whereas the 
net effect of very minor mergers and smooth gas accretion, 
even where it induces instabilities, may be to "spin up" the 
disks). 

As discussed in § 14.3.31 there is also an interesting 
regime of parameter space, namely minor mergers with 
mass ratios ~ 1:20- 1:10 or so, in which the characteristic 
merger timescales and secular/internal evolution timescales 
are comparable. The secondary may be large enough to 
induce a significant bar/non-axisymmetric response, but the 
merger/dynamical friction time may be sufficiently long that 
the primary could respond almost as if in isolation for several 
orbital periods. In such a case it becomes less clear whether 
the merger or the secular response of the disk is ultimately the 
dominant driver of evolution (and the answer probably de- 
pends on e.g. the exact orbital parameters and stability prop- 
erties of the disk, and may be sensitive to feedback, the gas 
phase structure and pressure support, and detailed halo struc- 
ture). In any event, it is clear that these processes require 
study in a more complete cosmological context, and can con- 
tribute significantly to the bulge population (especially in less 
bulge-dominated galaxies, below the typical thresholds we 
simulate) over a Hubble time of evolution. However, although 
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the bar itself may be triggered in the merger, the nature of the 
relative strength of the interaction and characteristic timescale 
for angular momentum loss make it not a violent process asso- 
ciated with the merger itself, but rather a secular process that 
should be considered more analogous to bars in non-merging 
systems. 

5. APPLICATION TO SEMI- ANALYTIC MODELS 

Our results clearly have potential uses as prescriptions for 
analytic and semi-analytic models of galaxy formation. Here, 
we summarize and give some simple recommendations for 
these applications. 

When a merger is identified in a semi-analytic model, the 
two key quantities we can predict here are the mass fraction of 
the disk that is destroyed (violently relaxed into a bulge) and 
the fraction of the cold gas in the disk that will lose angular 
momentum and contribute to the bulge by forming a compact 
starburst. 

First, the stellar disks: in a merger of secondary mass M2 
with primary mass M\ , the secondary is destroyed (adding M2 
to the bulge) and the mass within a radius enclosing w M2 in 
the primary is violently relaxed. If the primary were pure disk, 
this would add 2M2 to the bulge. However, one can imagine 
the limit where the primary is entirely bulge-dominated in- 
side that radius (with the stellar disk dominant only at much 
larger radii) - then the violent relaxation of the merger will 
act primarily to heat existing bulge stars, and only a mass 
1 M2 will be added to the bulge. Obviously, its also true that 
if the total disk mass of M\ is less than M%, then that is a 
maximum to how much can be added to the bulge (i.e. really 
MIN(M2, fdiskMi) is added). For most purposes, this factor 
2 possible range is not critical in the semi-analytic models, 
and picking a constant (effective mean) fraction (0-1) x M2 
to violently relax in the primary in all mergers is acceptable. 
However, if more detail is desired, an estimate of the mass 
profiles of bulge plus disk components in the primary can be 
used to determine the total primary disk mass within a ra- 
dius enclosing a mass « M2 , and then that w ill be the frac- 
tion violently relaxed. For a iHernquisll ( 1 19901) bulge and ex- 
ponential dis k obeying r oughl y the observed size-mass re- 
lations from IShen et all (12003b . the primary disk mass that 
should be violently relaxed in a merger with mass M2 can be 
approximated as fdi s k,*M 2 /(l + [Mi/M 2 ]"), where /disk.* = 
(1 ~/gas)(l - /bulge) is the mass fraction of the stellar disk (rel- 
ative to the baryonic galaxy) and the term (1 + [M\/M2\ a ) is a 
correction for e.g. the relative sizes of the two components as 
a function of mass ratio and other properties (for the assump- 
tions above, a w 0.3-0.6, depending on the details of the disk 
mass profile). 

Two clarifications should be emphasized. First, these 
derivations only apply to cases where the secondary is suf- 
ficiently massive that it survives to merge with the center of 
the primary. If the secondary is destroyed or shredded by tidal 
forces before merger, then it will not add either its own mass 
or any violently relaxed mass to the bulge. This generally 
occurs in the limits of smaller mass ratios (< 1 : 10, which 
we have considered), but is included in some models. Sec- 
ond, for most applications, the masses Mi and M2 should be 
taken to be the baryonic masses within the galaxies (stars in 
the galaxy and cold gas - not diffuse stellar halo or pressure- 
supported hot gas in the extended halo). This is how we have 
defined our models and fits to our simulations (although those 
simulations do include dark matter and extended halo gas and 
stars) throughout. The halos are much more extended, and 



much lower density, so they merge and mix more efficiently, 
and do not strongly participate in the central violent relaxation 
process that defines the bulge. Moreover, there can be a wide 
range in halo masses for galaxies of similar mass - but most of 
these halos are large and often independent substructures that 
should not be used to define e.g. the mass ratios of merging 
encounters. What dark matter is carried in with the galaxies is 
that enclosed in their stellar effective radii R e , which tends to 
track the baryonic mass much more closely than, say, the total 
halo mass, so it is not a bad proxy to still define mass ratios, 
etc. in terms of the baryonic masses. 

Next, in such an encounter, our analysis provides a means 
to estimate the fraction of the cold gas mass in the pre-merger 
stellar disks that should lose angular momentum and be fun- 
neled into a nuclear starburst. The cold gas inside some ra- 
dius Rgas/Rd will participate in this starburst, where R gsis is 
given by Equation (JVj. There are five variables that go into 
this equation: (1) / gas , which we define as the mass fraction of 
the disk that is in cold (rotationally supported) gas (i.e. if the 
disk is 50% cold gas, then regardless of the bulge fraction of 
the galaxy, / gas = 0.5. Note that we only care about cold, rota- 
tionally supported gas. Hot gas in the galactic halos can cool, 
of course, and form new stars, but that process is relatively 
independent of the merger, and is not related to angular mo- 
mentum loss (also because the hot gas is pressure-supported, 
it is fairly resistance to significant redistribution in the merger, 
and if anything will tend to be shocked to even higher temper- 
atures rather than forming stars in the short-lived merger). (2) 
/disk = (1 -/bulge), the total (gas plus stellar) baryonic mass 
fraction of the disk. (3) fj, = M2/M1, the mass ratio of the 
merger (defined as above). (4) 9 and b, equivalently the orbital 
parameters of the merger. As discussed in § 14.3.21 for cases 
that will merge the appropriate limit is £> — > 0, since most of 
the action will occur on the final merging passages after the 
angular momentum is removed. We discuss what should be 
adopted for the orbital inclination 9 below. (5) Rd, the scale 
length of the disk stars. 

In any semi-analytic or analytic model, variables (l)-(3) 
should be well-known beforehand. Given some choice of or- 
bital parameters and an assumed mass distribution of the disk, 
it is trivial then to translate Equation into a fraction of the 
gas that will burst. Because orbital parameters are generally 
undetermined in these models, there are two choices for the 
assumed orbital inclination 9. First, one could draw a random 
value of 9 for each merger (uniformly sampling in cos(#) as 
appropriate for an isotropic orbit distribution), and use Equa- 
tions <f9b-(TT0b for each merger. Alternatively, we can average 
over a random distribution of orbits and quote an "effective" 
orbital dependence F(9,b) for Equation ((9). Note that this 
is only strictly appropriate if all disks have the same mass 
profiles and those are such that the enclosed mass is linear 
in R/Rd (otherwise the appropriate average would have to be 
weighted by other terms such as (1 — / gas ) in Equation |7J. In 
any case doing so yields an effective mean orbital dependence 
F(0,fc)«1.2. 

The only remaining issue is the assumed mass profile of 
the disk. Here, models have some freedom. As we have em- 
phasized, the exact profile (e.g. choice of exponential disk or 
some other profile) does not have a dramatic effect. What is 
important, however, is the assumption of how the gas is dis- 
tributed relative to the stars. Recall, Equation Ql, with the 
variables above inserted, gives that the gas inside some ra- 
dius R g!ls = xRd (where x is a constant depending on those 
variables, and Rd is a characteristic scale length of the stellar 
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disk) should lose angular momentum and participate in the 
burst. Given a gas mass profile M gas {R/R e , g:ls ), in terms of a 
characteristic gas disk scale length ^ e ,gas, this gives the gas 
mass that bursts, M g:ls (xRa/R e . gas ). For our simulations, we 
have generally assumed (and can see that it is a good approx- 
imation) that the gas and stellar disks initially trace one an- 
other (7?d ~ R e , gas)- However, since our derivation and Equa- 
tion (0 show that it is the gas inside some fraction of the stel- 
lar disk half-mass radius R^ that loses angular momentum, 
then if the gas is e.g. much more extended than the stars, a 
lower gas fraction will end up in the burst. We discuss this 
in § 14.3.51 and consider how such situations may in fact arise 
owing to e.g. supernova feedback blowing gas out to large 
radii. Semi-analytic models therefore have some freedom in 
adopting these prescriptions based on their implicit assump- 
tions about feedback and disk formation, encapsulated effec- 
tively in our prescriptions as the ratio of the stellar to gas disk 
scale lengths Rd/R,,.^. Lacking some detailed model for both 
values in the semi-analytic models, a constant value ~ 1 is 
probably a good choice (with the exact choice reflecting im- 
plicit assumptions about feedback and outer disk formation). 

Those prescriptions define both the violently relaxed and 
starburst components induced in mergers of arbitrary mass 
ratios, gas content, and orbital parameters. If desired, appro- 
priate scatter (a factor ~ 2) can be added to both quantities, 
reflecting the scatter we see between various numerical real- 
izations (although it should still be ensured that, with scatter, 
the implied violently relaxed and burst fractions are within the 
sensible physical limits). 

Altho ugh not discussed here, in Hopki nset al.l 
(2008a.g.he.d) we consider how the sizes and velocity 
dispersions of these components should scale, and we refer to 
those papers for detailed analysis of those results. Briefly, we 
note that in the absence of dissipation, it is straightforward 
to calculate the size of the dissipationless component (the 
violently relaxed stars from the pre-merger stellar disk), 
given phase space and energy conservation. Roughly, this 
implies that the component will have the same (modulo 
projection effects since it transforms from a disk to a sphere) 
scale radius as the disk (or radius within the disk) from 
which it forms. Again, conservation of energy in subsequent 
dissipationless re-mergers, along with the assumption of 
pr eserved profile shape ( which we demonstrate is reasonable 
in iHopkins et al.l 1200 8 g) yields the evolution in subsequent 
events of these radii (in a re-merger of masses Mi and M2, 
the dissipationless bulge component will have final size 
R//R\ rs (1 + /i) 2 /(l + ^Rx/Ri)). Dissipation complicates 
this - it is possible to solve separately for the size of the 
dissipational component by allowing for energy loss in the 
collision followed by (after angular momentum loss) collapse 
to a self-gravitating limit, and then subsequently evolve the 
component as a dissipationless body, added with the violently 
relaxed co mponents to give a tot al bulge effective radius. 
Fortunately, Covington et al. (2008) perform such an exercise 
and we show in Hopkins et al.1 (1200 8 a) that their results can 
be conveniently approximated (in both an analytic manner 
and as a fit to the results of numerical simulations) by the 
scaling: fl e (bulge) = fl e (/ sb = 0)/(l +/ sb //o), where / sb is 
the total mass fraction of the bulge/spheroid which origi- 
nally formed dissipationally (as opposed to being violently 
relaxed), R e (f s b = 0) is the radius the system would have if 
purely dissipationless (calculated as described above), and 
fo ~ 0.25-0.35 is a constant. 

Our modeling could also be applied in the manner described 



in § l4.3.4l to fly-by (non-merging) encounters, but we caution 
that these are usually ill-de fined in semi-analytic models (and 
if adopted, the cautions in § !4.3.4l about the appropriate mean- 
ing of the impact parameter adopted should be borne in mind). 
In any case, the rapid suppression of bursts with increasing 
impact parameter means that such cases should be relatively 
unimportant in a representative cosmological ensemble. 

6. DISCUSSION AND CONCLUSIONS 

We have derived a general physical model for how disks 
survive and/or are destroyed in mergers and interactions. Our 
model describes both the dissipational and dissipationless 
components of the merger, and allows us to predict, for a 
given arbitrary encounter, the stellar and gas content of the 
system that will be dissipationlessly violently relaxed, dissi- 
pationally lose angular momentum and form a compact cen- 
tral starburst, or survive (without significant angular momen- 
tum loss or violent relaxation) to re-form a disk. We show 
that, in an immediate (short-term) sense, the amount of stel- 
lar or gaseous disk that survives or re-forms following a given 
interaction can be understood purely in terms of simple, well- 
understood gravitational physics. Knowing these physics, our 
model allows us to accurately predict the behavior in full hy- 
drodynamic numerical simulations across as a function of the 
merger mass ratio, orbital parameters, pre-merger cold gas 
fraction, and mass distribution of the gas and stars, in sim- 
ulations which span a wide range of parameter space in these 
properties as well as prescriptions for gas physics, stellar and 
AGN feedback, halo and initial disk structural properties, red- 
shift, and absolute galaxy masses. 

The fact that we can understand the complex, nonlinear be- 
havior in mergers with this analytic model, and moreover that 
(for given conditions at the time of merger) our results are in- 
dependent of the details of prescriptions for gas physics, star 
formation, and feedback, owes to the fact that the processes 
that strip angular momentum from gas disks and violently re- 
lax stellar disks are fundamentally dynamical. 

Gas, in mergers, primarily loses angular momentum to in- 
ternal gravitational torques (from the stars in the same disk) 
owing to asymmetries in the galaxy induced by the merger 
(on the close passages and final coalescence of the secondary, 
during which phase the potential also rapidly changes, scatter- 
ing and violently relaxing the central stellar populations of the 
stellar disk). 5 Hydrodynamic torques and the direct torquing 
of the secondary are second-order effects, and inefficient for 
all but pathological orbits. 

Once gas is efficiently drained of angular momentum, there 
is little alternative but for it to fall to the center of the galaxy 
and form stars, regardless of the details of the prescriptions 
for star formation and feedback - we show that even strong 
supernova-driven winds (with mass loading efficiencies sev- 
eral times the star formation rate and wind mass-loading ve- 
locities well above the halo escape velocity) do not signifi- 
cantly effect our conclusions. Such processes, after all, can 
blow out some of the gas, but they cannot fundamentally al- 
ter the fact that cold gas with no angular momentum will 
be largely unable to form any sort of disk, or the fact that a 
galaxy's worth of gas compressed to high densities and small 

5 We note again that although we have described these asymmetries as 
"bars" or "bar-like" at certain points in this paper, there are a number of 
prop erties of the non-ax isymmetric distortions induced in mergers (discussed 
in § 14.3.21 and § 14.3.71 that make them - at least over the short relaxation 
timescale of the merger - dynamically distinct from traditional bar instabili- 
ties in isolated systems. 
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radii will inevitably form a large mass in stars. 

For these reasons, many processes and details that are im- 
portant cosmologically (systematically changing e.g. the pre- 
merger disk gas fractions) - in some sense setting the initial 
conditions for our idealized study of what happens in mergers 

- do not alter the basic dynamical behavior within the mergers 
themselves, and therefore do not change our conclusions. 

Figure [16] summarizes our results for the ensemble of our 
simulations. We compare the fraction of the baryonic galaxy 
mass in the merger remnant that is in a surviving post-merger 
disk to that predicted by our simple model scalings, and find 
good agreement over the entire range in disk and bulge mass 
fractions sampled, with surprisingly small scatter given the 
complexity of behavior in mergers. We highlight several of 
the parameter studies, showing that - for fixed mass ratio, 
orbital parameters, and gas content at the time of the final 
merger, none of these choices systematically affect our pre- 
dictions (note that these are not the only parameters varied 

- the complete list is discussed in § [2] but it is representa- 
tive). That is not to say they cannot affect them indirectly, by 
e.g. altering how much gas is available at the time of merger 

- but it emphasizes that the processes we model and use to 
form our predictions, the processes that dominate violent re- 
laxation and the loss of angular momentum in gas in mergers, 
are fundamentally dynamical. 

This allows us to make robust, accurate physical predic- 
tions independent of the (considerable) uncertainty in feed- 
back physics and sub-resolution physics of the ISM. Regard- 
less of how those physics alter the "initial" conditions, they 
do not change basic dynamical processes, and so do not intro- 
duce significant uncertainties in our model. 

In turn, this means that we can use our model to under- 
stand just why and how feedback is important for the cosmo- 
logical s urvival of disks. Why, in short, have various works 
(see e.g. ISpringel & Hern quist 200H iRobertson et al.i r2006a; 
iGovernato et alJl2007h concluded that strong feedback is es- 
sential for enabling disk survival in mergers? Our results 
show that it is not that feedback somehow makes the disk 
more robust to the dynamical torques within the merger, in 
any instantaneous sense. These torques, at least within the 
critical radii where the gravitational perturbation from the 
merger is large and in resonance, are sufficiently strong that 
any reasonable feedback prescription is a dynamically neg- 
ligible restoring force. Rather, feedback has two important 
effects that fundamentally alter the conditions in the merger: 
first, it allows the galaxy to retain much higher gas content 
going into the merger. Without feedback from e.g. star for- 
mation and supernovae contributing to heating and pressur- 
izing the ISM and redistributing gas spatially, isolated gas- 
rich disks may be unstable to fragmentation. Even if frag- 
mentation is avoided, it is well-known that star formation in 
simulations proceeds efficiently under these conditions. This 
would leave the disks essentially pure stars (even for ideal- 
ized simulations beg inning with ~ 100% gas disks; see e.g. 
Springel et al. 2005b) by the time of the merger, which guar- 
antees that a major merger will inevitably violently relax the 
stars (this is a simple collisionless mixing process, and under 
such circumstances is inescapable). With large gas fractions, 
however, the system relies on stripping angular momentum 
from the gas to form new bulge stars, which in turn relies on 
internal torques from induced asymmetries in the stellar disk. 
If the gas fractions are sufficiently large, there is little stel- 
lar disk to do any such torquing, and the gas survives largely 
intact. 



Second, feedback from supernovae and stellar winds moves 
the gas to large radii, where it does not feel significant torques 
from the merger. Again, recall that the most efficient torquing 
is driven by the internal stellar disk of the galaxy, and as such 
is most efficient at torquing gas within small radii (this can be 
thought of as analogous to the well-known co-rotation condi- 
tion for isolated disk bars). If star formation-driven feedback 
has blown much of the gas to large radii, then there is little 
gas inside the radius where torques can efficiently strip an- 
gular momentum, yielding little induced starburst and largely 
preserving the gas disk at large radii. 

Not only can we qualitatively identify these requirements 
for feedback processes, but we can more precisely use our 
model to set quantitative limits on how much gas must be 
retained and/or the radii it must be redistributed to in order 
to enable disk survival under various conditions. This also 
clearly implies that disks must be able to avoid fragmentation 
and strong local gravitational instabilities when they achieve 
these gas fractions. This provides a valuable constraint for 
feedback models - how those models affect star formation ef- 
ficiencies, the "blowout" of gas, and the local hydrodynamic 
state (effective equations of state and phase structure) of ISM 
gas - and should be useful for calibrating their (still largely 
phenomenological) implementations in both numerical and 
semi-analytic models of galaxy formation. 

Our predictions are also of interest in any cosmological 
model for the emergence of the Hubble sequence, since they 
apply not just to disk-dominated galaxies but to small disks in 
bulge-dominated systems. We give a number of simple pre- 
scriptions for application of our conclusions to analytic and 
semi-analytic models of galaxy formation, which can be used 
to predict the distribution of bulge to disk ratios in cosmo- 
logical ensembles. But even without reference to a full such 
model, a number of interesting consequences are immediately 
apparent. 

First, it is a well-known problem that theoretical models 
systematically overpredict the abundance and mass fractions 
of bulges in (especially) low-mass galaxies. This is true even 
in e.g. semi-analytic models, which are not bound by resolu- 
tion requirements and can adopt a variety of prescriptions for 
behavior in mergers. However, it is also well-established ob- 
servationally that disk gas fractions tend to be very high in 
this regime, with l arge populations of g a s-dominated disks 
at < 10 10 Afo (iBell & de Jong! 1200 ll Ifannappanl 12004 
iMcGaughl 120051) . Our models predict that bulge forma- 
tion should, therefore, be strongly suppressed in precisely 
the regime required by observations. For e.g. disks with 
< 10 9 M© where observations suggest typical gas frac- 
tions ~ 60-80%, our results show that even a 1:1 major 
merger would typically yield a remnant with only ~ 30% 
bulge by mass - let alone a more typical 1:3-1:4 mass-ratio 
merger, which should yield a remnant with < 20% bulge. That 
is not to say that it is impossible to form a bulge-dominated 
system at these masses, but it should be much more difficult 
than at high masses, requiring either unusually gas-poor sys- 
tems, violent merger histories, or rarer merging orbits that are 
more efficient at destroying disks. Our conclusions therefore 
have dramatic implications for the abundance of bulges and 
typical morphologies and bulge-to-disk ratios at low galaxy 
masses and in gas-rich systems. Low-mass systems, when a 
proper dynamical model of bulge formation in mergers is con- 
sidered, should have lower bulge-to-disk ratios - by factors of 
several, at least - than have been assumed and modeled in pre- 
vious theoretical models. Whether this alone is sufficient to 
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FIG. 16. — Summary of our comparison between simulations and analytic model for the mass of disks in merger remnants as a function of appropriate orbital 
parameters, merger mass ratio, and pre-merger cold gas content. We plot our model prediction versus the simulation remnant disk fraction for all ~ 400 full 
hydrodynamic merger simulations considered in this paper (shown in both a linear and logarithmic scale). Symbols encode some of the parameter studies we 
consider: orbital parameters, galaxy masses, initial merger redshift, choice of feedback prescription, merger mass ratio, and presence or absence of black holes, 
as labeled. For each subset of simulations, we sample a wide range in initial and pre-merger gas fractions / gas = — 1 . Solid line is a one-to-one relation. In all 
cases, our predictions agree well with the simulations, with no systematic offsets owing to any of the parameters we have varied. At high /disk, our predictions 
are accurate to an absolute uncertainty ~ 0.05 - 0. 10 in /disk- At low /disk ^ 0. 1 , our predictions are accurate to a factor ~ 2 — 3 (down to /d; s k < 1%, where it is 
difficult to reliably identify disks in the remnant). 



resolve the discrepancies with the observations remains to be 
seen, but it is clearly of fundamental importance that future 
generations of models incorporate this scaling. 

Second, the importance of this suppression owing to gas 
content in disks will be even m ore significant a t high redshifts. 
Observations suggest (see e.g. Erb et al. 2006b) that by z ~ 2, 
even systems with masses near ~ L* (M* ~ 10 °— 10 u Mq) 
may have gas fractions as high as / gas ~ 0.6. In this regime, 
the same argument as above should apply, dramatically sup- 
pressing the ability of mergers to destroy disks. Moreover, 
since most of the mass density is near L*, this can change not 
just the behavior in a specific mass regime but significantly 
suppress the global mass density of spheroids, modifying the 
predicted redshift history of bulge formation. (Note that this 
will not change when stars form by very much, so it has little 
or no effect on e.g. the ages of z = spheroids). 

This redshift evolution may also explain the solution to a 
fundamental problem in reconciling observed disk popula- 
tions with CDM cosmologies. Integrated far enough back in 
time, every galaxy is expected to have experienced a signifi- 
cant amount of major merging. In extreme cases, the mass of 



the system when it had its last such merger may be so small 
that it would not be noticed today, but in general, it does not 
require going far back in redshift (to perhaps z ~ 2-4 before 
almost every z = galaxy should have had such a merger). 
How, then, can the abundance of systems with relatively small 
(or even no) visible bulges be explained? Our conclusions 
here highlight at least part of the answer: as you go back 
in time, the gas fractions of systems are also higher, nearing 
unity. So even though, integrating sufficiently far in time, ev- 
ery system has experienced major mergers, it is also true that 
the systems were increasingly gas-rich, and therefore that the 
impact of those mergers was more and more suppressed. Only 
mergers at later times, below certain gas fraction thresholds, 
will typically destroy disks. 

Third, to the extent that bulge formation is suppressed at 
incre asing redshifts, the exis tence of an M^- M\, u \ ge relation 
(e.g. iMagorrian et alj 119981) implies that black hole growth 
should also be suppressed. Indeed, bulge formation is sup- 
pressed specifically because gas cannot efficiently lose angu- 
lar momentum in mergers if the systems are gas-dominated - 
if the gas cannot lose angular momentum efficiently, then it 
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certainly cannot efficiently be accreted by the nuclear black 
hole. Since this pertains to gas on the scales of galactic disks, 
it is probably not relevant for the formation of "seed" black 
holes at very high redshift, but it will in general inhibit the 
growth of black holes owing to early merging activity. At the 
same time, of course, higher gas fractions in general imply in- 
creasing fuel supplies for black hole growth, so the effects are 
not entirely clear, and more detailed models are needed to see 
how this impacts the history of black hole growth and quasar 
luminosity functions. Nevertheless, this may in part explain 
why, above z ~ 2 (where, for the argument above, these effects 
become important for the global mass density of spheroids), 
the global rate of black hole growth (i.e. total quasar lumi- 
nosity density) appears to decline much more rapidly with in- 
cre asing redshift than the star formation rate density (compa re 
e.g. lHopkins & Beacomll2006l:lHopkins et alJl2008Al2007dl) . 

Fourth, our models imply that a large fraction of bulges and 
disks survive mergers together, rather than being formed en- 
tirely separately. It is often assumed that classical bulges - 
being similar to small ellipticals in most of their properties 
- were formed initially in major mergers, as entirely bulge- 
dominated systems, and then accreted new gaseous and stellar 
disks at later times. Although nothing in our modeling would 
prevent this from happening, our analytic and simulation re- 
sults generically lead to the expectation that a large (perhaps 
even dominant) fraction of the bulge population did not form 
in this manner, but rather formed in situ from minor mergers 
or less efficient major mergers (in e.g. very gas-rich systems). 
Observations tracing the evolution of disk components, kine- 
matics, and morphology in the last ~ lOGyr increasingly sug- 
gest that such co -formation or disk regeneration scenario is 
common (see e. g. Hamm er et alj|2005t IConselice et alj|2005t 
iFlores et"aT1l2006t iPuech et alJl2008T and references therein). 
In short, a system with a mass fraction ~ 0.1-0.2 in a bulge 
could be the remnant of an early, violent major merger (when 
the system was ~ 0.1 times its present mass) with a re- 
accreted disk, or could be the remnant of a typical (low to 
intermediate gas fraction) 1:10-1:5 mass ratio minor merger, 
or could even be the remnant of a gas-rich major merger (mass 
ratio < 1 : 3, if / gas is sufficiently large). 

Based on a simple comparison of typical merger histories, 
we would actually expect that the minor merger mechanism 
should be most common, but all may be non-negligible. 
Fundamentally, the physics forming the bulge (torquing the 
gas within some radius owing to internal asymmetries and 
violently relaxing stars within a corresponding radius) are 
the same in all three cases, and moreover other indicators 
such as their stellar populations will be quite similar (in all 
cases, the bulge will appear old: this is both because the 
central stars in even present-day disks are much older than 
those at more typical radii, and because in any case star 
formation will cease within the bulge itself, as opposed to 
the ongoing star formation in the disk, and stellar population 
age estimates are primarily sensitive to the amount of recen t 
or ongoing star formation; see e.g. Trag er et all [2000). 
This is also not to say that mergers are the only means of 
producing bulges. Secular evolution of e.g. barred disks 
probably represents an increasingly important channel for 
bulge evol ution in later-type and mo r e gas-rich systems 
(see e.g. IChristodo ulou et al l Il995ah ISheth etalJ 120031: 
Mayer & Wadslevl I2004T Debattist aet alj 12004 Uogee et all 
2004t iKormendv & Kennicuttl 12004 iMarinova & Joged 
2007), and may even be related (albeit through longer 
timescales of "isolated," post-merger evolution and different 



physics) to initial bar formation or "triggering" in mergers. 
More detailed theoretical work and analysis of cosmological 
simulations is needed to develop observational probes that 
can distinguish between these histories. 

Further work is specifically needed to investigate the pro- 
cesses at work in minor mergers with mass ratios ^1:10 
(fi ~ 0.05-0.1), which cosmological simulations suggest are 
an important contrib utor to the growth of disks, especially in 
later-type systems jMaller et aLH2006l: iFakhouri & Mall200l 
IStewart et al.ll2.Q07h . In more minor mergers /x <C 0.1, the sec- 
ondaries are sufficiently small and dynamical friction times 
sufficiently long that the disk is unlikely to feel significant ex- 
ternal perturbations. More major mergers fj, > 0.1, the cases 
of interest here, induce sufficiently large responses in the disk 
and evolve sufficiently rapidly that they can be considered 
"merger-dominated" for the reasons in § 14.3.21 & 14.3.71 But 
in the intermediate regime, internal amplification of instabili- 
ties in a traditional secular fashion may occur on a timescale 
comparable to or shorter than the evolution of the secondary 
orbit, potentially leading to a more complex interplay between 
the two. It is not entirely clear whether such a system would 
remain "locked" to the driven perturbation, or function as a 
purely secular system (merely initially driven by the presence 
of the secondary), or some nonlinear combination of both. A 
more detailed comparison of the relevant timescales for these 
processes and their relation to e.g. cosmological triggering of 
bars and large-scale non-axisymmetric modes in disks will be 
the subject of future study (in preparation). 

Our results are also of direct interest to models of spheroid 
formation in ellipticals and SO galaxies. As discussed in §Q] 
it is increasingly clear that embedded sub-components - con- 
stituting surviving gaseous and stellar disks - are both ubiq- 
uitously observed and critical for theoretical models to match 
the de tailed kinematics a nd isophotal shap es of observed sys- 
tems (iNaab et al.1 120061: ICox et al.1 [2006aHbl: iRobertson et alj 
2006b; Uesseit et alJl2007b iHopkins et al.ll2008allbh . We have 
developed a model that allows us to make specific predic- 
tions for how disks survive mergers, including both the sur- 
vival of some amount of the pre-merger stellar disks and the 
post-merger re-formation of disks and rotationally supported 
components from gas that survives the merger without losing 
most of its angular momentum. 

Figure [16] shows that we can extend these predictions with 
reasonable accuracy to surviving rotational systems contain- 
ing as little as ~ 1 % of the remnant stellar mass, comparable 
to small central subcomponents and subtle feat ures giving rise 
to e.g . slightly disky isophotal shapes (see e.g.lFerrarese et al.1 
11994 lLauer et alJ 120051 iMcDermid etail 120061) . Owing to 
the combination of resolution requirements and desire to un- 
derstand the fundamental physics involved, most theoretical 
studies of these detailed properties of ellipticals have been 
limited to idealized studies of individual mergers. Our re- 
sults allow these to be placed in a more global context of cos- 
mological models and merger histories. Moreover, our mod- 
els allow the existence of such features (or lack thereof) to 
be translated into robust constraints on the possible merger 
histori es and gas-richness o f spheroid-forming mergers. Fur- 
ther, [Hopkins et al. ( 2008a, b), studied how the dissipational 
starburst components arising in gas-rich mergers are critical 
to explaining the observed properties and scaling relations of 
ellipticals, and how these components can both be extracted 
from and related to observed elliptical surface brightness pro- 
files. Because both the starburst and surviving disks arise 
from gas in mergers, the combination of constraints from the 
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central stellar populations, studied therein, with constraints on 
the survival and/or loss of gas angular momentum in mergers 
studied here, should be able to break some of the degenera- 
cies in e.g. pre-merger gas fractions and merger histories in 
order to enable new constraints and understanding of spheroid 
merger histories, and new tests of models for spheroid forma- 
tion in gas-rich mergers. 

These points relate to a number of potentially testable pre- 
dictions of our models. These include the in situ formation of 
bulges from various types of mergers, and possible associated 
stellar population signatures, the presence of embedded disks 
in ellipticals, and how their sizes and mass fractions scale with 
e.g. the masses and formation times of ellipticals (and how 
this relates to gas fractions and stellar populations in observed 
disks). In general, for similar merger histories, the increas- 
ing prevalence of later type galaxies (SO's and SOa's) at lower 
masses where disks are characteristically more gas rich is a 
natural consequence of our predictions here, and it is straight- 
forward to convert our predicted scalings into detailed pre- 
dictions for the abundance and mass fractions of disks given 
some simplified merger histories. To the extent that these pro- 
cesses also give rise to disk heating and/or increasing veloc- 
ity dispersions in disks, or changing kinematics in both disks 
and bulges, then there should be corresponding relationships 
between galaxy shapes, kinematics, and bulge-to-disk ratios 
along the Hubble sequence. We investigate these possible cor- 
relations and tests in subsequent papers (in preparation). 

Altogether, our results here elucidate the relevant physics 
important for both dissipational and dissipationless bulge for- 
mation in mergers. They support a new paradigm in which to 
view bulge and disk formation: gas-richness is not simply a 
"tweak" to existing models of bulge formation and disk de- 



Abadi, M. G., Navarro, J. R, Steinmetz, M., & Eke, V. R. 2003, ApJ, 597, 21 
Athanassoula, E. 2002a, ApJ, 569, L83 
— . 2002b, Ap&SS, 281, 39 
— . 2003, MNRAS, 341, 1179 
— . 2005, MNRAS, 358, 1477 

Athanassoula, E., Bienayme, O., Martinet, L., & Pfenniger, D. 1983, A&A, 
127, 349 

Athanassoula, E., Lambert, J. C, & Dehnen, W. 2005, MNRAS, 363, 496 
Athanassoula, E., & Misiriotis, A. 2002, MNRAS, 330, 35 
Balcells, M., Graham, A. W., & Peletier, R. F. 2007, ApJ, 665, 1084 
Barazza, F. D., Jogee, S., & Marinova, I. 2008, ApJ, 675, 1194 
Barnes, J. E. 1988, ApJ, 331, 699 

Barnes, J. E. 1998, in Saas-Fee Advanced Course 26: Galaxies: Interactions 
and Induced Star Formation, ed. R. C. Kennicutt, Jr., F. Schweizer, J. E. 
Barnes, D. Friedli, L. Martinet, & D. Pfenniger, 275-+ 

Barnes, J. E., & Hernquist, L. 1996, ApJ, 471, 115 

Barnes, J. E., & Hernquist, L. E. 1991, ApJ, 370, L65 

Barton, E. J., Arnold, J. A., Zentner, A. R., Bullock, J. S., & Wechsler, R. H. 

2007, ApJ, 671, 1538 
Bell, E. E, & de Jong, R. S. 2001, ApJ, 550, 212 

Bell, E. F, Phleps, S., Somerville, R. S., Wolf, C, Borch, A., & 

Meisenheimer, K. 2006, ApJ, 652, 270 
Bender, R., Burstein, D., & Faber, S. M. 1992, ApJ, 399, 462 
Bender, R., Surma, P., Doebereiner, S., Moellenhoff, C, & Madejsky, R. 

1989, A&A, 217, 35 
Benson, A. J. 2005, MNRAS, 358, 551 

Berentzen, I., Athanassoula, E., Heller, C. H., & Fricke, K. J. 2003, MNRAS, 

341, 343 
— . 2004, MNRAS, 347, 220 

Berentzen, I., Shlosman, I., Martinez- Valpuesta, I., & Heller, C. H. 2007, 
ApJ, 666, 189 

Binney, J., & Tremaine, S. 1987, Galactic dynamics (Princeton, NJ, Princeton 

University Press, 1987) 
Borderies, N., Goldreich, P., & Tremaine, S. 1989, Icarus, 80, 344 
Bournaud, E, & Combes, F. 2002, A&A, 392, 83 
Bournaud, F, Jog, C. J., & Combes, F. 2005, A&A, 437, 69 



struction in mergers. Rather, if disks are sufficiently gas rich, 
the qualitative character of mergers is different, with ineffi- 
cient angular momentum loss giving rise to disk-dominated 
remnants. This process is not inherently governed by poorly- 
understood feedback physics (although such feedback may be 
critical for establishing the conditions necessary in the first 
place), but rather by well-understood gravitational physics, 
and as such is robust and fundamentally inescapable. As- 
pects of galaxy populations such as the continuum of relative 
bulge and disk mass ratios are not simply consequences of 
e.g. different amounts of accretion, but can arise owing to the 
continuum in efficiencies of disk destruction as a function of 
merger mass ratios, orbital parameters, and gas content. The 
relative (lack of) abundance of bulges at low galaxy masses 
and high redshift is a basic consequence of the dynamics of 
how gas loses angular momentum in mergers, even for simi- 
lar merger histories. In short, the baryonic physics of mergers 
ensures that, despite the near self-similarity of the physics and 
merger histories of their host halos, disk and bulge formation 
are not a self-similar process, influenced dramatically (well 
out of proportion to the absolute cold gas mass fractions) by 
the gas-richness of the baryonic systems. 
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